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Introdiiccion 



Los chorros Hquidos son usados en diversos procesos industriales y maquinarias cuando 
se desea, entre otros objetivos, aumentar la relacion superficie - volumen, dirigir el liquido 
hacia un bianco determinado, transferir energia entre el liquido y un solido, o la formacion de 
un hilo a partir de un material fundido. 

En algunas de estas aplicaciones se desea que el chorro rompa en gotas, mientras que en 
otras se busca que el flujo permanezca en su configuracion inicial. La presencia en el flujo de 
perturbaciones libres o forzadas (por diversos medios), que crecen en el espacio y en el 
tiempo, es la causa de la inestabilidad del mismo y de su consecuente ruptura. El estudio de 
inestabilidades constituye la herramienta que permite estimar si la ruptura se produce o no, 
dados los parametros de interes correspondientes al escurrimiento, y en el caso de producirse 
la ruptura predecir las caracteristicas de las gotas que se forman. 

El caracter inestable de un flujo puede ser clasificado de diversas maneras de acuerdo a 
la evolucion que muestran las perturbaciones. Una de estas posibles clasificaciones es la de 
flujo convectiva o absolutamente inestable, en donde para el primer caso la perturbacion crece 
alejandose aguas abajo de la fuente de perturbacion, y en el segundo crece aguas abajo y 
aguas arriba de la fuente. 

El analisis de la evolucion de una perturbacion infinitesimal en el espacio y en el 
tiempo, es segiin algunos autores el medio adecuado para establecer la anterior clasificacion, y 
es para este caso en donde se aprecia su relevancia. De tal estudio surgen predicciones tales 
como la velocidad de propagacion y los limites del paquete de ondas resultante de la 
descomposicion de la perturbacion, y la tasa de crecimiento espacio - temporal. 

La excitacion de la superficie del chorro por medio de fuerzas electricas, proporciona un 
medio de control para la evolucion de las perturbaciones, pudiendose gobemar el caracter 
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inestable del mismo, los parametros de las gotas formadas luego de la ruptura, y 

eventualmente la trayectoria de las mismas. 

Para el presente trabajo, de particular importancia son las modificaciones que por este 
medio de excitacion, pueden producirse para el caso de un chorro convectivamente inestable, 
sometido a perturbaciones tipo pulso excitadas acusticamente, en terminos de la tasa de 
crecimiento y de la expansion de los limites del paquete de ondas. Con este trabajo se 
pretende realizar un analisis experimental de la dinamica de este tipo de escurrimientos. Los 
resultados de este estudio, revisten importancia tanto desde un punto de vista practico como 
fundamental. 

Existe un modelo de transicion a la turbulencia en capas limites, que posee un 
tratamiento matematico similar al que aqui se trata. La verificacion experimental de ese 
modelo es sumamente dificultosa por las caracteristicas inherentes al escurrimiento, hecho 
que obliga a emplear tecnicas experimentales muy complejas. En los chorros liquidos los 
estudios experimentales son relativamente mas sencillos, por lo tanto a partir de sus resultados 
podrian inferirse algunas caracteristicas de las inestabilidades en el escurrimiento antes 
mencionado. 

Desde el punto de vista de las aplicaciones de la excitacion electrica a los chorros 
liquidos, se puede considerar la posibilidad de controlar a voluntad la transicion entre las 
inestabilidades convectivas y absolutas. Este fenomeno guarda una estrecha relacion con el 
cambio de regimen entre un chorro establecido en su configuracion inicial y en regimen de 
goteo. Si bien aqui se analiza un chorro convectivamente inestable, los resultados que se 
espera obtener para el comportamiento de los limites del paquete de ondas frente a la 
excitacion electrica, estan fuertemente relacionados con el fenomeno aludido. En 
consecuencia, con este trabajo se pretende conseguir resultados que permitan avanzar en el 



8 



Introdiiccion 



desarroUo de un dispositivo de control que posibilite el pasaje a voluntad de un regimen a 
otro. 

El presente informe se ha organizado de la siguiente manera: 

En el capitulo 1 se definen los chorros Hquidos y se trata someramente sobre su 
estabilidad, caracterizando los regimenes, modos y rupturas en gotas que pueden presentarse 
cuando son excitados en forma libre o forzada. Las vias de excitacion forzada son descriptas, 
finalizando con el tratamiento de algunas aplicaciones practicas. 

La descripcion de un trabajo teorico desarroUado con anterioridad a este, y que es su 
marco de referenda, es el objeto del capitulo 2. Previo a ello, el concepto de inestabilidad y su 
estudio es revisado y ampliado, describiendo las distintas clasificaciones que pueden darse. El 
capitulo finaliza con la presentacion en forma concisa del objetivo de esta tesis. 

En el capitulo 3 se describe el principio de funcionamiento y los diferentes componentes 
del montaje experimental, que ha sido desarroUado y construido en ocasion de este trabajo. 

La presentacion de los resultados experimentales, y la comparacion con las previsiones 
teoricas, es efectuada en el capitulo 4. Los limites del estudio, y el tratamiento de las senales 
para la determinacion de los parametros de interes, son previamente introducidos. 

Esta tesis concluye con la exposicion de sus resultados mas notorios, junto a algunas 
consideraciones que tal vez podrian ser tenidas en cuenta en estudios posteriores. 
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1.1 - EVTRODUCCION 

Este capitulo tiene como objeto presentar el sujeto del trabajo, los chorros Uquidos. Para 
ello se comienza por dar una definicion de chorros y de su estabilidad, haciendo luego una 
descripcion de chorros no excitados que incluye una resena historica de trabajos de 
investigacion sobre este campo. Posteriormente se realiza una descripcion de chorros 
excitados aciistica y electricamente. Se mencionan finalmente las principales aplicaciones 
practicas de los chorros dando una breve explicacion de cada una. 

1.2 - DEFINICION DE CHORRO 

Un chorro es un escurrimiento resultante de la inmersion de un fluido que surge desde 
un orificio, dentro de otro fluido. Para este trabajo, el fluido que surge del orificio sera un 
Hquido, el orificio circular, y el otro fluido sera el aire sujeto a presion atmosferica. EI 
escurrimiento sobre el cual se trata aqui esta dentro de la categoria de flujo de superficies 
libres, que se caracterizan por la puesta en movimiento de medios continues en los cuales 
coexisten varias fases liquidas o gaseosas separadas unas de las otras por interfases. 

1.2.1 - Estabilidad de un chorro 

Las leyes de conservacion de magnitudes mecanicas y termodinamicas tales como la 
masa, la cantidad de movimiento, o la energia, permiten establecer un sistema de ecuaciones 
diferenciales que describen la dinamica del flujo. 

La resolucion analitica de este sistema de ecuaciones generales de movimiento para un 
fluido en regimen estacionario, con las condiciones de borde correspondientes, solo puede 
lograrse en contados casos. Este hecho obhga a trabajar en general con un sistema de 
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ecuaciones mas sencillo obtenido mediante la introduccion de hipotesis simplificativas, para 
hallar asi una solucion aproximada de la dinamica del escurrimiento. Es evidente que la 
introduccion de las hipotesis mencionadas hace que la solucion aproximada sea valida dentro 
de ciertos Hmites impuestos per ellas. Sin embargo, aiin si se conoce la solucion exacta del 
sistema de ecuaciones original, no significa que el movimiento descripto por la teoria 
realmente exista en la naturaleza. De hecho la presencia de perturbaciones al escurrimiento no 
tenidas en cuenta en el modelo, alteran la dinamica del flujo y pueden conducir a su 
desestabilizacion. 

Un chorro esta compuesto de una parte aproximadamente cilindrica y de otra parte en la 
cual aparecen gotas que pueden formar una nube, segiin su numero y tamano. La parte 
cilindrica del chorro corresponde en general a la solucion aproximada (por ejemplo a un perfil 
constante de velocidad en un cilindro infinito). En la segunda parte, las diferencias que 
pueden ser observadas entre la solucion aproximada y la experiencia (ruptura del chorro y 
aparicion de gotas), indican claramente el caracter inestable del chorro que no logra 
amortiguar el crecimiento de las perturbaciones que conducen a la ruptura. 

Un estudio de la estabilidad de un chorro consiste entonces en primer termino de la 
obtencion de un perfil de velocidades sobre un dominio con fronteras fijas que es la 
configuracion inicial de equilibrio, y en segundo termino del analisis de la evolucion de dicho 
perfil cuando se impone una perturbacion sobre la frontera. Un fiujo sera inestable si existe al 
menos una perturbacion cuya amplitud aumenta con el tiempo. 

La Uamada ruptura de un chorro es un fenomeno de desintegracion en gotas que se da 
al cabo de una cierta longitud del escurrimiento. Las caracteristicas de este fenomeno 
dependen del diametro del chorro, de su velocidad, y de las propiedades fisicas del fluido del 
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chorro y la atmosfera circundante. Los diferentes tipos de ruptura de un chorro seran 
descriptos mas adelante. 

1.2.2 - Chorros sin excitacion y con excitacion 

En los chorros sin excitacion, las perturbaciones iniciales no son impuestas 
extemamente. Las fluctuaciones en la geometria del chorro cuando el mismo se forma en el 
orificio aparecen asociadas a efectos naturales aleatorios como variaciones radiales del flujo, 
turbulencia, vibraciones, vortices, etc., que genericamente son calificados como ruido. 

Los chorros liquidos excitados son aquellos en los cuales las perturbaciones iniciales a 
la forma de la superficie del chorro son preferenciales, es decir promovidas extemamente. La 
excitacion preferencial de un chorro liquido puede ser efectuada por medio acustico, electrico, 
y tambien combinaciones de ambos. 

1.3 - CHORROS LlOUroOS SIN EXCITACION 

En este apartado se detallan los aspectos salientes de este tipo de chorros, comenzando 
por los tipos de ruptura que pueden presentarse de acuerdo a la velocidad del escurrimiento 
dando lugar a diferentes regimenes, presentando luego las curvas de estabilidad, para luego 
hacer una resena historica de trabajos de investigacion sobre inestabilidades en chorros no 
excitados, y finalmente describir los distintos modos de perturbacion de superficie de chorros. 

1.3.1 - Tipos de ruptura en funcion de la velocidad 

Un chorro no se desintegra en gotas en forma inmediata a la salida. En general, entre el 
orificio de salida del chorro y la parte en donde comienza la aparicion de gotas, existe una 
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zona cuya longitud es Uamada longitud de chorro intacto, caracterizada por poseer una 
superficie aproximadamente cilmdrica, alterada por la presencia de perturbaciones. Reitz y 
Braco [1], establecieron una clasificacion de los chorros basandose en la longitud intacta del 
chorro y en la relacion entre el diametro del orificio de salida y el diametro medio de las 
gotas, distinguiendo los cuatro regimenes siguientes, que pueden verse en la figura 1.1. 

Regimen de Rayleigh: en este regimen, el diametro de las gotas, para una ruptura 
normal, es mayor que el diametro del orificio. Este es un regimen que aparece para 
velocidades bajas. Las fuerzas ligadas a la tension superficial en el chorro son quienes 
favorecen el aumento de la amplitud de las perturbaciones en el tiempo y provocan la ruptura 
del chorro. En este tipo de chorros, la suma de las superficies de las gotas producidas es 
inferior a la superficie del cilindro liquido del cual provienen. Este es el regimen que tiene 
interes a los efectos del presente trabajo. 

Primer Regimen Aerodinamico: la velocidad en este caso es mas importante que en el 
anterior. Se observa un niimero importante de gotas de menor diametro que el del orificio. 
Este regimen se caracteriza por el crecimiento de perturbaciones de longitudes de onda 
mayores que el radio del orificio. Aqui las fuerzas asociadas al movimiento relativo entre el 
chorro y la atmosfera circundante, y aquellas debidas a la tension superficial, son las que 
provocan el crecimiento de las perturbaciones. 

Segundo Regimen Aerodinamico: en este regimen, las gotas formadas son aiin de 
tamano mas chico que en el caso anterior, y la longitud intacta del chorro es importante. La 
velocidad del chorro es considerable, y el regimen es normalmente turbulento. La formacion 
de pequenas gotas es el resultado del crecimiento de inestabilidades de longitudes de onda 
menores sobre la superficie del chorro. En este regimen, la suma de las superficies de las 
gotas es mayor a la del cilindro liquido del cual provienen. Las fuerzas que surgen de la 
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interaccion del chorro con el medio que lo circunda son principalmente las responsables de la 
desestabilizacion del chorro, mientras que las fuerzas asociadas a la tension superficial se 
oponen al crecimiento de las perturbaciones. 
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Regimen de Rayleigh 




0 



Primer Regimen Aerodinamico 




Regimen de Atomizacion 



Figura 1.1 —Diferentes regimenes de ruptura de un chorro. 
Regimen de Atomizacion: los tamanos de las gotas formadas son aun menores que los 
regimenes precedentes, y la longitud intacta del chorro es practicamente nula. La velocidad 
del chorro es mucho mayor. En este regimen, las gotas resultantes de las perturbaciones de 
longitud de onda menores son arrancadas de la superficie del chorro formando un cono 
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alrededor del chorro que se desintegra. Las perturbaciones que tienen una longitud de onda 
grande pueden sin embargo progresar hasta el final del chorro resultando en gotas de diametro 
importante. 

Ranz [2], propuso otra clasificacion de los regimenes de chorros, basandose en que la 
fisica de los dos regimenes aerodinamicos anteriores presenta suficientes similitudes como 
para reunirlos en uno solo. La clasificacion la establecio en ftincion de los numeros de 
Ohnesorge (tambien Uamado "numero de estabilidad", Oh) y de Reynolds {Re), defmidos por: 

Oh= 

Ml 

donde /ui es la viscosidad dinamica del liquido, pi su densidad, a es el diametro del 
chorro, /la tension superficial entre el liquido y el gas, y Uo es la velocidad media del chorro. 




100 1000 10000 100000 

Re 

Figura 1.2 - Diferentes tipos de ruptura de un chorro en funcion de los numeros de 

Ohnesorge (Oh) y de Reynolds (Re). 
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1.3.2 - Longitud de ruptura de un chorro liquido 

Puede definirse la longitud de ruptura de un chorro como aquella en la cual el chorro es 
continuo. Esta longitud depende del regimen del escurrimiento. Numerosos investigadores 
han representado las curvas de variacion de la longitud de ruptura en funcion de la velocidad 
del chorro (Uamadas curvas de estabilidad). Aqui se presenta una curva tipica en la figura 1.3. 

Bower et al [3], para las velocidades mas bajas, observaron una primera region de 
escurrimiento en donde toma lugar el regimen de goteo, y a partir del punto A, un chorro 
establecido. La zona entre los puntos A y B es la correspondiente al regimen de Rayleigh, en 
la cual la longitud de ruptura aumenta linealmente con la velocidad. El primer regimen 
aerodinamico se situa entre los puntos B y C, y esta caracterizado por una disminucion de la 
longitud de ruptura cuando la velocidad del chorro aumenta. El segundo regimen 
aerodinamico esta comprendido entre los puntos C y D observandose un aumento de la 
longitud de ruptura con la velocidad. A partir del punto D, la longitud de ruptura decrece 
hasta un minimo con el aumento de la velocidad, en correspondencia con el regimen de 
atomizacion. En esta ultima zona, la determinacion de la longitud de ruptura se hace dificil ya 
que los metodos opticos no pueden ser empleados y las sondas que pueden utilizarse tocan el 
chorro y perturban asi la medida. 

Para completar la descripcion anterior, Arai et al [4], determinaron que si la presion 
ambiente es superior a la atmosferica normal, las longitudes de ruptura se ven disminuidas. 
Ademas, si la velocidad es grande, la longitud de ruptura tiene dependencia de la forma del 
inyector. Finalmente, si el fenomeno de cavitacion se presenta dentro del inyector, se observa 
en general una disminucion de la longitud de ruptura, resultado obtenido por Hiroyasu et al 
[5]. 
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Longitud de ruptura 




Atomizacion 



Rayleigh 



Velocidad del chorro 



Figura 1.3 - Esquema de una curva de estabilidad de un chorro. 



1.3.3 - Resena historica del estudio de inestabilidades en chorros sin excitacion 



A continuacion se presenta en forma resumida la evolucion de los primeros trabajos de 
investigacion sobre chorros liquidos no excitados, para lo cual se ha prestado atencion a la 
resena que figura en el trabajo de McCarthy y MoUoy [6]. 

Probablemente los trabajos pioneros sobre la descripcion del fenomeno de ruptura de 
chorros son los que fueron publicados por Bidone (1829) [7] y por Savart (1833) [8]. El 
primero de ellos estudio la geometria de chorros formados en orificios de secciones 
transversales regulares diferentes de la cilindrica. Posteriores investigaciones sobre este 
aspecto fueron hechas por Magnus (1855) [9], Buff (1856) [10] y Boussinesq (1877) [11]. 
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Savart proporciono los primeros datos cuantitativos sobre la ruptura de chorros, 
proponiendo dos leyes experimentales, a saber: para un diametro constante, la longitud de 
ruptura de un chorro es directamente proporcional a la raiz cuadrada de la presion que 
provoca el escurrimiento (es decir proporcional a la velocidad); y para velocidad constante, la 
longitud de ruptura es proporcional a su diametro. Savart presento sus resultados 
experimentales empleando las Uamadas curvas de estabilidad, y fue el primero en utilizar tales 
curvas. 

El primer trabajo teorico sobre estabilidad de chorros fue efectuado por Plateau (1856) 
[12], quien mostro que una columna cilindrica estatica de liquido puede ser inestable si su 
longitud excede su perimetro ya que, si esto fuera asi, podrian formarse dos gotas cuya area 
superficial total seria menor que la del cilindro a partir del cual se formaron, y entonces esta 
ultima configuracion tendria una energia superficial menor. El trabajo de Plateau explico en 
parte los resultados de Savart, y fue la base para el desarroUo de la teoria de las 
inestabilidades en chorros. 

1.3.3.1 - Chorros no viscosos en el vac/o 

Rayleigh (1878) [13], estudio la inestabilidad de chorros no viscosos en el vacio, 
situacion que se aproxima, por ejemplo, a un chorro de agua de baja velocidad en aire, en 
donde el escurrimiento esta sujeto a fuerzas de inercia y a aquellas asociadas a la tension 
superficial. Basandose en la suposicion de que el chorro estaba perturbado infinitesimalmente, 
encontro que el chorro solamente es inestable cuando las perturbaciones son axisimetricas (es 
decir, unicamente fiincion de la coordenada axial del chorro), y si las longitudes de onda de 
las perturbaciones son mayores que el producto 2m, donde a es el radio del chorro en el 
estado no perturbado. Asumiendo que la amplitud de la perturbacion crece en forma 
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exponencial con el tiempo, determino la tasa de crecimiento en funcion de la longitud de 
onda, su valor maximo y el valor de longitud de onda para el cual se presenta. Ademas logro 
predecir la longitud de ruptura del chorro suponiendo que la misma se produce cuando la 
amplitud de la perturbacion de mayor tasa de crecimiento alcanza el valor del radio del chorro 
sin perturbar. De esta manera determino que la longitud de ruptura varia proporcionalmente a 
la velocidad del chorro y a la relacion ln(a/5o), con So la amplitud inicial de la perturbacion. 
Las conclusiones a las que arribo Rayleigh surgen de un analisis lineal, ya que al desarroUar la 
forma de la superficie del chorro perturbado en series de Fourier se quedo solamente con los 
terminos de primer orden. Por lo tanto, los resultados de esta teoria solo son aplicables cuando 
las amplitudes de las perturbaciones son mucho menores que el radio del chorro no 
perturbado. Obviamente esta limitacion no flie tenida en cuenta en la derivacion de la 
expresion para la longitud de ruptura, ya que en la zona de ruptura la amplitud de la 
perturbacion tiene el mismo valor que el radio del chorro. No obstante, en primera instancia, 
un analisis lineal constituye una buena aproximacion ya que la amplitud de las perturbaciones 
toma importancia solo en el 10 o 20 % final de la longitud del chorro antes de su ruptura. 

Trabajos posteriores de muchos investigadores confirmaron la teoria de Rayleigh entre 
ciertos limites. Uno de ellos es el de Smith y Moss (1916-1917) [14], quienes usando curvas 
de estabilidad mostraron que la longitud de ruptura en un rango de velocidades moderadas es 
proporcional a la velocidad del chorro, destacando tambien que la relacion ln(a/5o) (la escala 
relativa de la perturbacion inicial) permanecia constante con un valor de 13 para todos los 
liquidos empleados en sus estudios experimentales. Para velocidades por encima del rango 
mencionado, observaron una disminucion de la longitud de ruptura ante el aumento de la 
velocidad. Ademas observaron que el primer maximo de la curva de estabilidad 
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(correspondiente al limite superior de velocidad del rango de proporcionalidad) esta 
emparentado con la aparicion de turbulencia en el escurrimiento. 

1.3.3.2 - Efecto de la viscosidad 

Rayleigh (1879) [15] tambien investigo acerca de chorros viscosos sometidos a la 
accion de la tension superficial. La ecuacion para la tasa de crecimiento resultante de su teoria 
no es provechosa dada su complejidad. Esta teoria indica que cuando las fuerzas viscosas son 
mucho mayores que las fuerzas de inercia, la longitud de onda para la maxima tasa de 
crecimiento tiende a infinito, por lo cual puede decirse que un efecto del incremento de la 
viscosidad es incrementar la longitud de onda para la cual la tasa de crecimiento es maxima. 

Un posterior analisis de la estabilidad de chorros viscosos en el vacio fue hecho por 
Weber (1931) [16], obteniendo resultados similares a los de Rayleigh pero mas sencillos 
debido a la introduccion de algunas simplificaciones en el tratamiento matematico. De la 
teoria de Weber surgen expresiones que permiten conocer la maxima tasa de crecimiento y la 
longitud de onda correspondiente a ella, ademas de la longitud de ruptura del chorro en 
funcion de las fuerzas de inercia, las viscosas y las de tension superficial. El efecto de la 
viscosidad, de acuerdo a este analisis, es el de amortiguar las perturbaciones haciendo 
disminuir la tasa de crecimiento efectiva. 

Merrington y Richardson (1947) [17], experimentalmente comprobaron los resultados 
de Weber en cuanto a la linealidad entre la longitud de ruptura y la velocidad del chorro en un 

rango de velocidades determinado, a partir del cual la relacion se invierte pasando a disminuir 
la longitud de ruptura con el aumento de la velocidad, y encontrando que el limite superior del 
rango de proporcionalidad crece junto con el aumento de la viscosidad. Meister y Scheele 



21 



Generalidades sohre chorros 



(1969) [18] determinaron que la teoria de Weber se circunscribe a valores del numero de 

Ohnesorge que cumplen con la siguiente inecuacion: 

— <0,2 
Oh 

En la mayoria de las situaciones esta inecuacion requiere que la viscosidad del chorro 
sea mayor que 160 cP, aspecto que fue ignorado aparentemente por muchos investigadores, 
Weber incluido. 

1.3.3.3 - Efecto de la velocidad del chorro y de la interaccion con el medio circundante 

Una explicacion del decrecimiento de la longitud de ruptura con el aumento de la 
velocidad fue propuesta por primera vez en el trabajo de Eisenklam y Hooper (1958) [19]. 
EUos encontraron que la relajacion del perfil de velocidad de un chorro laminar de alta 
velocidad (perfil de velocidad parabolico), hacia un perfil de velocidad constante (cuyo valor 
es el mismo que la velocidad media del perfil parabolico), luego de la formacion del chorro en 
el orificio, trae aparejada una redistribucion de energia lo suficientemente violenta como para 
que el chorro se desintegre. La conclusion a que conduce el trabajo es que los chorros con 
perfiles parabolicos de velocidad en su inicio son mas inestables que aquellos en donde el 
perfil de velocidad es tipo "tapon". Este tipo de estudios fiie retomado posteriormente por 
Leib y Goldstein (1986) [20], quienes concluyeron que los chorros con perfil de velocidad 
parabolico son mas estables (tanto sea en terminos de tasa de crecimiento como de rango de 
frecuencias inestables), que aquellos en donde el perfil de velocidad es tipo tapon. Sin 
embargo, las perturbaciones que son neutrales cuando el perfil de velocidad es parabolico 
pueden finalmente volverse crecientes en la medida que el perfil de velocidad se relaja hacia 
uno del tipo tapon. 
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En el campo de los chorros en regimen turbulento, uno de los primeros trabajos 
corresponde a Grant y Middleman (1966) [21]. EUos mostraron que cuando el chorro es 
turbulento, la longitud de ruptura vuelve a crecer con el aumento de la velocidad, y 
propusieron una correlacion para ella en funcion de la relacion entre las fuerzas de inercia y 
las de tension superficial, que no es lineal. 

En terminos del tamano de gotas producidas, el primer trabajo fue efectuado por 
Scheubel (1927) [22]. Posteriormente Schweitzer (1937) [23], Uego a la conclusion de que 
las teorias que justifican la atomizacion de chorros de alta velocidad por la accion del aire 
circundante solamente, son invalidas. El considero que la turbulencia del escurrimiento, 
inducida tal vez en el orificio, juega un papel preponderante en el proceso de atomizacion, 
aunque fuertemente ayudada por la friccion con el aire. Estudios posteriores entre los cuales 
se encuentra el de Ranz (1958) [2], retomaron la idea de que la fi-iccion con el aire es la causa 
de la ruptura en gotas para chorros en este regimen. 

El primer trabajo teorico acerca de chorros liquidos en un medio liquido circundante, 
fue realizado por Tomotika (1935) [24], quien empleo un analisis lineal para desarroUar una 
relacion entre la tasa de crecimiento de perturbaciones axisimetricas y la longitud de onda 
para sistemas no miscibles en general. Sus resultados son aplicables cuando la relacion entre 
la viscosidad de los liquidos toma un valor finito. Es interesante destacar que si la viscosidad 
del liquido y del medio se desprecia, se Uega a la misma expresion de Rayleigh para chorros 
no viscosos en el vacio, y si la viscosidad del chorro no puede ser ignorada pero si la del 
medio circundante, se Uega a la misma expresion de Weber. La teoria de Tomotika ha sido 
confirmada en parte y ampliada por muchos trabajos posteriores, entre los cuales se 
encuentran los de Miesse (1955) [25], Ranz y Dreier (1964) [26] y Meister y Scheele (1969) 
[18]. 
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1.3.3.4 - Efectos no lineales 

Las teorias lineales predicen que los chorros sujetos a la accion de la tension superficial 
se desintegran en gotas de tamano uniforme, con una gota formada por longitud de onda de 
perturbacion, y que la forma de la superficie del chorro perturbada corresponde a una 
sinusoide que se amplifica. Las observaciones experimentales sin embargo muestran que 
sobre todo cerca de la zona de ruptura, la forma del chorro es diferente a la sinusoidal, y 
destacan la presencia de gotas pequenas entre gotas de mayor tamano. Para explicar estos 
fenomenos se han desarrollado trabajos experimentales y teorias no lineales por parte de 
numerosos investigadores entre los que se encuentran Donnelly y Glaberson (1966) [27], 
Wang (1968) [28], y Yuen (1968) [29]. 

A titulo de resumen de los estudios de inestabilidades en chorros sin excitacion, puede 
decirse que una gran parte de ellos se basan en general en una linealizacion del problema y 
que las predicciones que se efectuan surgen de analizar las curvas caracteristicas que vinculan 
el numero de onda k con la frecuencia oj. De estas curvas, la perturbacion que resulta como 
dominante es la que se corresponde con la maxima tasa de crecimiento. Esta perturbacion se 
toma como referente para la determinacion de longitudes de ruptura y tamanos de gota. Los 
efectos de los diferentes parametros fisicos sobre estas magnitudes han sido largamente 
estudiados y las corroboraciones experimentales han sido en general satisfactorias. 

1.3.4 - Descripcion de los diferentes modos que pueden ocurrir en chorros sin excitacion 

Se ha mencionado en los parrafos precedentes que la teoria lineal de la propagacion de 
perturbaciones en chorros es bien conocida y ha sido tratada por numerosos autores [30, 31, 
32, 33]. La situacion mas general que ha sido estudiada se corresponde a un chorro con 
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perturbaciones axiales y angulares. La superficie perturbada puede describirse 
matematicamente como sigue, para el caso mas general: 

r(z, e,t) = a + Tio- (jj) 

En donde a es el radio del chorro sin perturbar, ijo^s la amplitud inicial de la 
perturbacion, z es la coordenada axial, 0 la coordenada angular, t el tiempo, n es un entero 
igual o mayor que cero, k es el denominado numero de onda, que vale: k = 2-n / X , con X la 
longitud de onda de la perturbacion axial, y <» la frecuencia. 

En el caso de los chorros sin excitacion, el numero de onda es real, y la frecuencia es un 
numero complejo. Esta situacion sera analizada luego en el capitulo siguiente. 

La figura 1.4 muestra los modos mas comunes que pueden desarroUarse en un chorro 
circular. El primer modo corresponde a n=0 y k?^, en el cual el chorro es perturbado en forma 
axisimetrica, que fue denominado "varicoso" por Rayleigh, en alusion a la similitud con una 
vena varicosa, y puede ser encontrado ademas en la literatura con los nombres "salchicha" y 
"simetrico". La superficie del chorro perturbada con este modo se caracteriza por sucesivas 
contracciones y expansiones. Para el presente trabajo, este es el modo que reviste interes, 
hecho que sera explicado posteriormente. 

El segundo modo es aquel en el cual n=l, y en este caso el desplazamiento de la 
superficie del chorro respecto de la posicion de equilibrio es proporcional al coseno de la 
coordenada angular. La seccion transversal del chorro, mientras permanece circular, es 
desplazada en relacion al eje del chorro, adquiriendo forma de onda. Este modo se denomina 
comunmente "sinuoso". 
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Los modos siguientes, para n=>2, son conocidos genericamente con el nombre de 
"flauta". Para n=2, la seccion transversal del chorro es eliptica, y los diametros mayor y 
menor cambian de posicion entre ellos altemativamente. 




Figura 1.4 - Esquema de varios modos de perturbaciones en la superficie de chorros. [34] 

1.4 - CHORROS LIQUIDOS EXCITADOS ACUSTICAMENTE 

En este caso las perturbaciones iniciales son impuestas al escurrimiento empleando 
actuadores piezoelectricos, que son forzados a vibrar en frecuencias que comprenden el rango 
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acustico y eventualmente el ultrasonico. Como resultado del movimiento vibratorio del 
piezoelectrico, se inducen oscilaciones en la salida del chorro, que promueven la aparicion de 
las perturbaciones sobre el chorro. El actuador piezoelectrico puede ser montado de diferentes 
maneras, segiin sea que el chorro se forme con un inyector o un capilar. En el caso de 
formarse el chorro con un capilar, el actuador piezoelectrico se coloca en general solidario a 
el. Cuando se emplea un inyector el actuador se coloca en el extremo de una camara 
especialmente conformada en el inyector (la perturbacion de presion en este caso viaja hasta 
el orificio), o bien el actuador puede ser usado directamente como soporte para proveer el 
orificio en donde el chorro es formado. 

En esta tesis se emplea la excitacion acustica combinada con la electrica, y como un 
adelanto a la descripcion del dispositivo experimental puede decirse que el montaje del 
actuador piezoelectrico corresponde al segundo caso antes mencionado. 

1.4.1 - Ruptura normal 

En chorros de baja velocidad, la ruptura normal es el resultado de la amplificacion de 
perturbaciones axisimetricas que desestabilizan el escurrimiento. Como se ha afirmado 
anteriormente, para chorros no excitados, la longitud de onda de tales perturbaciones debe 
estar dentro de un rango determinado, y las causantes de la desintegracion son aquellas cuyas 
longitudes de onda estan cercanas a la que produce la maxima tasa de crecimiento de la 
amplitud de la perturbacion. La ruptura normal de chorros excitados o no, es similar, 
diferenciandose principalmente en que en el caso no excitado se observa una irregularidad en 
el tamano de gotas y en la separacion entre ellas. 

La ruptura normal en el caso de un chorro excitado esta gobemada por la frecuencia de 
excitacion impuesta al chorro, y no por la correspondiente a la maxima tasa de crecimiento 
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como en el caso de chorros no excitados. Se caracteriza por producir una gota por ciclo de la 
oscilacion, gotas de igual tamano, cuyo diametro resulta mayor que el correspondiente al 
chorro sin perturbar, y depende de la velocidad y radio del chorro, y de la longitud de onda de 
la perturbacion. 

La figura 1.5 muestra la ruptura normal de un chorro de agua desmineralizada formado 
con un capilar de 148 \\m de diametro intemo, fluyendo a una velocidad de 6,13 m/s, excitado 
con una senal sinusoidal de frecuencia 8036 Hz. La longitud intacta observada es 0,417 cm. 



Figural.5 - Ruptura normal de un chorro de agua desmineralizada. [35] 
1.4.2 - Ruptura con gota satelite (efectos no lineales) 

Este tipo de ruptura, que puede calificarse como anormal, se caracteriza por producir 
mas de una gota por ciclo de la oscilacion. Esta gota extra, a menudo no deseada, es 
generalmente mas pequena (aunque no siempre) que la gota principal, y son Uamadas 
satelites. 

La figura 1.6 muestra el comportamiento de un chorro liquido para diferentes numeros 
de onda adimensionales K = k-a , tomada de [36]. 
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Kf= 0-683 



Figura 1.6-Ruptura con gota satelite de un chorro de agua. 
Wetsel [37], reporta en su trabajo que la aparicion de gotas satelites se debe a efectos no 
lineales, y que este fenomeno esta ligado a una determinada banda de frecuencias de 
excitacion, a la escala relativa de la perturbacion inicial, al contenido de armonicos de la 
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frecuencia de excitacion en las perturbaciones y a la aparicion de ondas secundarias entre las 
crestas de las ondas de la perturbacion inicial. 

1.4.3 - Ruptura anormal (efectos no lineales) 

Otros tipos de ruptura anormal aparte de la anterior, estan descriptas en el trabajo de 
Bassani y Scaricabarozzi [35]. EUos distinguen los dos casos siguientes. 

1.4.3.1 - Ruptura con multiples sotas 

Este tipo de ruptura se presenta en las figuras a continuacion: 




Figura 1.8 -Ruptura anormal con multiples gotas de un chorro de agua desmineralizada. 
Velocidad: 4,01 m/s, frecuencia: 1873 Hz, longitud intacta: 0,53 cm. 
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Figura 1.7 -Esquema de ruptura anormal con multiples gotas. 
Cerca del punto de ruptura del chorro, se observan pequenas gotas por unidad de 
longitud de onda, que se unen para formar otras mas grandes. Esta region tiene una extension 
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igual a 2A, (Region i de la figura 1.7). A esta region le sigue otra de longitud variable entre una 
y nueve longitudes de onda, que presenta dos gotas por X (Region ii). Suficientemente lejos 
del extremo del capilar, se observan gotas equidistantes de igual diametro (Region iii). El 
tamano y la distancia entre gotas se ajusta al de la ruptura normal, en este ultimo caso. 

El analisis efectuado por los autores no explica la causa de la ruptura en pequenas gotas 
ni el mecanismo de su posterior coalescencia. Sin embargo, el analisis de la perdida de 
energia superficial propuesto, proporciona algunos elementos para predecir su aparicion. 

1.4.3.2 - Ruptura con multiples Glas de sotas 

Aqui luego de la ruptura se forman dos filas divergentes de gotas. La distancia entre 
gotas de la misma fila es igual a la de la ruptura normal, mientras que el diametro de las gotas 
resulta menor que en la ruptura normal y puede ser determinado analiticamente a partir de un 
balance de masas. El comienzo del proceso de ruptura es similar al caso anterior de ruptura 
con multiples gotas, en tanto que la divergencia del patron de gotas previene su coalescencia. 
Esta separacion se produce por la oscilacion transversal del tubo capilar. 




Figura 1.9 — Ruptura anormal con multiples filas de gotas de un chorro de agua 
desmineralizada. Velocidad: 4,01 m/s, fiecuencia: 545 Hz. 
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1.4.4 - Inestabilidades en chorros excitados acusticamente 

Desde un punto de vista de la teoria de inestabilidades, en este caso el problema es 
diferente con respecto al caso anterior. En tanto que antes se consideraba como referente a la 
perturbacion cuya longitud de onda se correspondia con la tasa de crecimiento maxima, aqui 
lo es aquella que tiene la frecuencia con que se realiza la excitacion. Este problema es 
conocido como "problema de la senal" (signalling problem). 

1.4.5 - Diferentes modes en chorros con excitacion acustica 

La superficie del chorro en este caso tambien puede ser descripta matematicamente por 
la ecuacion 1 . 1 dada para el caso de un chorro no excitado. La diferencia se establece en que 
para chorros excitados acusticamente la frecuencia es un numero real y el numero de onda es 
un numero complejo. Los modos que aparecen en este tipo de excitacion son entonces los 
mismos que pueden verse en la figura 1 .4. 

En sintesis, para chorros con excitacion acustica, la ruptura se produce para una 
perturbacion de longitud de onda que es fijada por la excitacion. La ruptura normal da como 
resultado gotas de tamano uniforme en cantidad de una por longitud de onda de perturbacion, 
y puede darse el caso de ruptura anormal cuando los fenomenos no lineales toman 
importancia, sobre todo en la zona final de la longitud de chorro intacto, a consecuencia de lo 
cual surgen gotas de diferente tamano o bien se forman filas divergentes de gotas. La 
situacion de ruptura anormal generalmente es no deseada en aquellos casos en donde se quiere 
controlar la trayectoria de gotas por medio de campos electricos, situacion que se presenta en 
algunas aplicaciones practicas, como se vera posteriormente. Sin embrago, en ocasiones la 
gota satelite es buscada para poder manipular una gota de menor tamano de la que surgiria de 
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una ruptura normal. El hecho de que las perturbaciones tienen una longitud de onda dada por 
la excitacion, establece diferencias en el estudio de inestabilidades frente a los chorros no 
excitados, en lo que hace al numero de onda y la frecuencia, pero no hace diferencias respecto 
de los modos de las perturbaciones. 

1.5 - CHORROS LIQUIDOS EXCITADOS ELECTRICAMENTE 

La existencia de una diferencia de potencial electrico entre un chorro liquido y aquello 
que lo rodea da origen a un campo electrico. Este campo electrico, dependiendo de las 
tensiones aplicadas, podra ser altemo o continuo. En este trabajo se prestara atencion al caso 
de campos continuos, no sin antes mencionar que la aplicacion de campos altemos puede 
devenir en la excitacion preferencial del chorro en forma similar a la excitacion aciistica. 

La existencia de un campo electrico continuo puede derivar en una modificacion de las 
caracteristicas del escurrimiento y de su estabilidad. La importancia de los efectos de tal 
excitacion es funcion de las caracteristicas del chorro (velocidad, conductividad electrica del 
liquido, etc.) y de los parametros del sistema de electrificacion empleado (longitud del 
electrodo, intensidad de campo, etc.) 

Las principales consecuencias posibles de la accion de un campo electrico continuo 
sobre un chorro son: la migracion de cargas hacia la superficie del chorro y la creacion de un 
chorro que se denomina electrificado, la modificacion de su estabilidad, y finalmente la 
excitacion del chorro con una perturbacion preferencial, diferente a la del caso de excitacion 
acustica. 

La migracion de cargas hacia la superficie del chorro, provocada por el campo electrico, 
es el resultado de la accion de las fuerzas coulombianas sobre los iones del liquido. Esta 
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migracion permite la obtencion de chorros electrificados que luego de su ruptura resultan en 
gotas cargadas. 

1.5.1 - Distintas formas de excitacion electrica 

Un chorro puede ser electrificado de diversas maneras. Bailey [38, 39] agrupo esos 
diferentes metodos de electrificacion de la siguiente forma. 

1.5.1.1 - Excitacion por conduccion 




Electrode 
conectado a tierra 



Liquido sometido a 
tension 



Figura 1.10 - Esquema de un sistema de excitacion por conduccion. 
Este primer sistema esta representado en la figura anterior. El polo positivo de una 
fiiente de alta tension esta unido al inyector en el interior del cual se encuentra el liquido. El 
otro polo esta en contacto con un electrodo hueco a traves del cual el liquido pasa sin tener 
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contacto. De esta forma el inyector se coloca a un potencial elevado y las cargas son 
inyectadas en el Hquido hasta que el potencial Uega a ser uniforme, por conduccion en el seno 
Hquido. La forma y la posicion del electrodo conectado a tierra permite controlar el campo 
electrico. En ciertos dispositivos el electrodo no existe, el campo electrico se establece entre 
el chorro y los cuerpos que lo rodean, por lo que ellos mismos se establecen directa o 
indirectamente como el potencial de tierra. 

1.5.1,2 - Excitacidn por inHuencia 



Liquido conectado 




Electrodo 
sometido a 
tension 



Figura 1.11 -Esquema de un sistema de excitacidn por influencia. 
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Este sistema se muestra en la flgura 1.11, y es el que reviste importancia para el 
presente trabajo. El sistema es analogo al precedente. El Hquido en este caso atraviesa un 
electrodo hueco sometido a un alto potencial sin tener contacto con el. Esta vez el chorro es 
conectado a tierra y el campo electrico creado entre el electrodo y el chorro atrae algunas 
cargas hacia la superficie del escurrimiento. El campo actiia aqui como separador de cargas 
presentes en el liquido. El liquido neutro de partida contiene, natural o artificialmente, un 
numero igual de cargas positivas y negativas que son repelidas o atraidas por el campo 
exterior. 

1.5.1,3 - Excitacion por inyeccion de carsa 



Punto de alta 
tension 



Zona de inyeccion 
de cargas 






Gotas cargadas 





Figura 1.12 - Esquema de un sistema de excitacion por inyeccion de carga. 
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Este es un sistema que se emplea cuando se trabaja con liquidos que son pobres 
conductores de la electricidad. Se electrifica el chorro por medio de la inyeccion de cargas 
con la ayuda de una punta ubicada en el interior del inyector y conectada al polo positivo de 
una fuente de alta tension, mientras que el inyector es conectado a tierra. Los iones creados 
son removidos por el chorro que escurre desde el inyector. 

1.5.1.4 - Excitacion por impacto de iones 




Gotas no cargadas 



Punto corona 



Gotas cargadas 



Figura 1.13 - Esquema de un sistema de excitacion por impacto de iones. 
Este ultimo sistema esta esquematizado en la figura 1.13. Consiste en el aporte de iones 
sobre la superficie del chorro o sobre las gotas al hacer pasar el liquido a traves de una region 
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en la cual se provoca una descarga corona. Los iones producidos por la descarga chocan con 
el Iiquido y son incorporados por el. 

1.5.2 - Inestabilidades en chorros excitados electricamente 

El estudio de inestabilidades en chorros fuertemente electrificados tiene una elevada 
diversidad de posibilidades y una descripcion de todas ellas escapa, en general, al objetivo de 
este trabajo. La cantidad de casos en que estos estudios han sido efectuados es ademas escasa, 
por lo cual se centrara la atencion sobre inestabilidades en chorros electrificados en ocasion 
de presentar los aspectos teoricos que hacen al marco del estudio experimental, dentro del 
capitulo siguiente. 

1.5.3 - Regimcncs cn chorros electrificados 

No existe un acuerdo generalizado acerca de la nomenclatura a seguir para clasificar los 
diferentes regimenes que pueden darse en chorros excitados electrostaticamente. A 
continuacion se expone una clasificacion que fue establecida en el trabajo de Cloupeau y 
Prunet-Foch [40]. Esta clasificacion no toma en cuenta aquellos casos en que el 
escurrimiento es forzado por medios neumaticos u otros, en los cuales la aplicacion del campo 
electrico tiene el objeto de cargar las gotas y de reducir ligeramente su diametro medio, sino 
aquellos en los cuales el campo electrico es la causa de la atomizacion del Iiquido en finas 
gotas. 

1.5.3.1 - Resimen de soteo 

En este regimen, la emision de gotas puede ocurrir a intervalos regulares de tiempo, sin 
la creacion de gotas satelite, de forma tal que las gotas tienen todas el mismo tamano. 
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La figura 1.14 representa esta situacion y muestra en la curva dos puntos, A y B, en los 
cuales se observa respectivamente cuando el menisco esta en la posicion mas elevada y 
cuando la gota se separa. En general las gotas son de mayor tamano que el del capilar, lo que 
Ueva a la emision de gotas grandes a bajas frecuencias. 



Figura 1.14 - Regimen de goteo sin formacion de satelites. Relacion entre la frecuencia de 

desprendimiento de gotas y la tension [40]. 

Contrariamente al caso anterior, en donde no hay formacion de satelites, cada gota 

puede estar acompanada por la formacion de uno o varios satelites. La figura 1.15 muestra 

que las gotas principales estan unidas por un largo y delgado filamento cuando se separan del 

capilar, el cual se rompe en un gran numero de pequenas gotas. 



Figura 1.15 - Regimen de goteo con gota satelite. Gotas principales unidas por un filamento. 
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Si el campo electrico en el extremo de la gota pendiente es suficiente para crear un 
chorro siguiendo el mismo proceso que para el modo cono - chorro descripto a continuacion, 
un gran numero de pequenas gotas son emitidas durante el ciclo de formacion de cada gota 
principal. 




O 



Figura 1.16 -Regimen de goteo con formacion de chorro. 
1.5.3.2 - Resimenes con formacion de conos de Taylor 

En estos casos el menisco toma la forma de un cono extendido en su extremo por un 
chorro permanente que rompe luego en gotas. Es de mucho interes porque pueden ser 
obtenidos para un amplio rango de conductividad electrica del liquido, pueden producirse 
aerosoles con un amplio rango de tamano medio de gotas incluido el rango submicronico, y 
los aerosoles son muchas veces monodispersos. 

1.5.3.2.1 - Regimen de cono-chorro 

La figura 1.17 muestra esquematicamente este regimen. Para liquidos de relativamente 
alta conductividad, la zona de formacion del chorro se limita al extremo del menisco, el cual 
puede tener generatriz recta (a) o curva (b). 

Para valores menores de conductividad, la zona de aceleracion se extiende mas hacia la 
base del cono. En el limite, esta zona puede comenzar en la salida del capilar. Sin embargo el 
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perfil del Hquido a la salida del capilar a menudo permanece similar a un cono, mas o menos 
abierto (c, d). 

Tambien puede darse una asimetria en el menisco (e). En el limite, el desprendimiento 
del chorro puede Uegar a estar en la periferia del extreme del capilar. 




b e d « 



Figura 1.17 — Diferentes formas del menisco en el regimen cono - chorro. 
Si la carga del chorro no es muy elevada, la ruptura del mismo en gotas toma lugar de la 
misma manera que los chorros neutros sometidos a inestabilidades varicosas, y aqui tambien 
es posible observar la aparicion de gotas satelite. Dependiendo de las condiciones, los 
aerosoles producidos pueden ser mono o poli dispersos. (Figura 1.18 a). 




> b 

Figura 1.18 - Ruptura del chorro en el regimen cono - chorro. Inestabilidades varicosas (a), 

inestabilidades sinuosas (b). 

Cuando la presion electrostatica en el chorro supera la presion capilar, las gotas mas 

grandes tienen una carga que supera la carga limite de Rayleigh, y entonces emiten un chorro 

que se desintegra en finas gotas. Si la tension es ligeramente aumentada, aparece la 
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inestabilidad sinuosa, y se incrementa a medida que el campo lo hace. El chorro rompe en este 
caso en finas gotas de muchos tamanos. (Figura 1.18 b). 

Algunas variantes del modo cono chorro se muestran en la siguiente figura. El caso (a) 
se presenta para valores de campo algo menores de los que establecen un chorro continuo, es 
un regimen de transicion. Los casos (b) y (c) se dan para campos electricos grandes. 




a b e 



Figura 1.19 - Variantes del regimen cono - chorro. Cono - chorro pulsado (a), 

multichorro (b). 

1.5.3.2.2 - Resimen de micro soteo 




Figura 1.20 - Regimen de micro goteo. 
En este caso la emision de gotas toma lugar de una en una como en el modo de goteo 
convencional, pero las gotas tienen un diametro mucho menor que el de la seccion de salida 
del capilar. En general para este regimen, en el extremo de un menisco conico o redondeado 
aparece un chorro de longitud variable, en el extremo del cual hay una acumulacion de liquido 
que concluye en la formacion de una gota. Luego de la separacion de la gota, el chorro toma 
una extension variable y el ciclo se reinicia. Este regimen aparece superpuesto al de cono - 
chorro para valores bajos de tension y frecuencia de goteo. 
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1.5.3.3 - Resimen de chorros simple y ramiffcado 

En ausencia de un campo electrico, el Hquido forma un chorro permanente si el caudal 
es lo suficientemente alto, por encima de un caudal critico. En este caso la aplicacion de un 
campo electrico de bajo valor tiene el efecto de cargar y acelerar el chorro ya establecido. 

El chorro tambien puede ser obtenido con caudales menores al critico, si la tension 
aplicada excede un cierto valor minimo. 

En cualquiera de los dos casos anteriores, cuando la tension es incrementada, el 
comportamiento del chorro al principio es similar al que se observa en el regimen cono - 
chorro, y luego se producen ramificaciones en niimero variable. 

La figura 1.21 muestra esquematicamente este regimen y sus diferentes aspectos. En el 
caso (a), el chorro rompe en grandes gotas posiblemente acompanadas de gota satelite. En (b), 
las gotas separadas emiten a veces uno o dos chorros delgados que rompen en finas gotas. En 
(c) se aprecia que la ruptura se produce por la amplificacion de inestabilidades sinuosas. En el 
caso (d) el chorro forma una sucesion de regiones engrosadas desde las cuales emergen uno o 
varios chorros delgados. En (e), el chorro se ensancha a poco de dejar el capilar y muchos 
chorros finos son emitidos alrededor. Finalmente en (f), para muy altos valores de tension, la 
parte permanente del chorro desaparece y el aspecto del aerosol cambia constantemente. 




Figura 1.21 - Chorros simple y ramificado. 
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1.5.3.4 - Resimen "ahusado" 

En los regimenes precedentes, la produccion de gotas toma lugar desde la ruptura de un 
chorro o directamente de la separacion de parte del liquido acumulado en el menisco. Un caso 
especial lo constituye aquel en el cual los dos tipos de formacion de gotas coexisten. 

Este regimen se representa en la figura 1.22, la cual muestra las diferentes etapas del 
ciclo de formacion de gotas. Un menisco redondeado (1) se alarga (2) y un chorro delgado 
emerge de su extremo rompiendo luego en gotas (3). La elongacion del menisco se acentua 
(4,5) y una estrangulacion ocurre en la base de la parte ahusada (6). Un volumen de liquido 
con forma de huso finalmente se separa (7). Luego se contrae formando una gota principal 
(8), mientras que el chorro remanente continua rompiendo en finas gotas (9,10). 
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Figura 1.22 - Representacion esquemdtica del regimen ahusado. 
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1.5.3.5 - Pasaje de uno a otro resimen 

En general, demasiados parametros interactuan (viscosidad del liquido, constante 
dielectrica, conductividad electrica, tension superficial, tipo de capilar, entre otros), y no 
siempre es posible conocer de antemano el regimen a obtener cuando se conocen los valores 
de los principales parametros. Asi, para un mismo sistema y un mismo liquido, pequenas 
diferencias en detalles experimentales como la forma y la capacidad de mojado del extremo 
del capilar, definen la iniciacion de un modo sobre otro y cambios en las caracteristicas del 
aerosol formado. 

En sintesis, diferentes regimenes tienen lugar en la excitacion electrica de chorros, 
siendo en algunos casos sutiles las diferencias que permiten clasificarlos, lo cual hace que 
puedan aparecer solapados cuando se toman en cuenta relaciones entre los parametros de 
interes. Las diferencias que se observan para cada modo en lo que respecta a los mecanismos 
de ruptura, promueven que las caracteristicas de las gotas formadas, como por ejemplo el 
tamano, scan diferentes para cada uno de ellos. 

1.6- INTERES DEL ESTUDIO SOBRE CHORROS. APLICACIONES 

Las aplicaciones practicas de los chorros Hquidos son numerosas y variadas, y este 
hecho, junto con el natural interes desde el punto de vista teorico, hace que el estudio de la 
inestabilidad de chorros adquiera gran relevancia. 

Se ha dejado en claro que la amplificacion de cierto tipo de perturbaciones conduce a la 
desestabilizacion de los chorros, lo que acarrea la posterior ruptura en gotas. Ligados a las 
aplicaciones practicas aparecen dos problemas bien diferenciados: por un lado puede interesar 
que el chorro no se desintegre en gotas (por ejemplo algunos procesos de Uenado de 
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recipientes y formacion de fibras metalicas), y por otro lado interesa que la ruptura se 
produzca dando como resultado gotas de un tamano determinado (procesos de atomizacion). 

La aplicacion de campos electricos a los chorros interesa desde el punto de vista de 
modificar la dinamica del chorro y de guiar las gotas cargadas que se producen a traves de una 
determinada trayectoria. 

A continuacion se describen de forma breve algunas de las aplicaciones de los chorros 
[30,38,39,41]. 

1.6.1 - Formacion de gotas en el proceso Sol-Gel 

Una aplicacion importante de los chorros liquidos es la desintegracion del chorro para 
formar gotas de liquido que contienen combustible nuclear. Estas gotas son luego secadas y 
sinterizadas para formar microesferas ceramicas para usar como combustible en reactores 
nucleares. Los mecanismos para la produccion de gotas son varios, pero siempre empleando 
chorros de baja velocidad [30]. 

Una solucion acuosa que contiene el metal a ser usado como combustible (uranio o 
plutonio) es dispersada en gotas dentro de un agente desecante liquido de naturaleza organica, 
el cual remueve el agua de la solucion formando una suspension coloidal semi solida, que es 
luego secada y sinterizada en homos. 

1.6.2 - Formacion de fibras metalicas 

Los chorros de metal liquido son usados para producir fibras metalicas que tienen 
aplicacion industrial o en joyeria. Una de las tecnicas mas simples es aquella en la cual el 
metal fundido es forzado a escurrir verticalmente a traves de un orificio, adquiriendo un 
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movimiento de caida libre, siendo solidificado durante el mismo. En este caso no se desea que 
el chorro rompa en gotas, y el calor de solidificacion debe ser removido antes de que se 
vuelva inestable. 

1.6.3 - Pulverizacion de pesticidas 

Esta aplicacion se presenta en la agricultura como medio de control de plagas. Los 
pesticidas son atomizados de diferentes maneras, con equipos portatiles manuales, equipos 
tirados por tractores o desde el aire empleando aviones. En forma convencional, la 
atomizacion se produce empleando diversos dispositivos, como por ejemplo pistolas 
atomizadoras, discos y copas rotantes, etc. En estos casos una gran parte del chorro se pierde a 
causa del movimiento de las gotas en el aire, dando como resultado una deposicion ineficaz 
del pesticida sobre el cultivo. Para aumentar la eficiencia se emplean dispositivos de 
atomizacion electrostatica, para los cuales esta demostrado experimentalmente que aumenta el 
rendimiento de la deposicion al generarse gotas cargadas electricamente. 

1.6.4 - Pulverizacion de pinturas 

La aplicacion de pinturas atomizadas permite mejorar el rendimiento del recubrimiento 
obteniendo en general espesores mas delgados y uniformes frente a otros metodos de 
aplicacion. Los pulverizadores convencionales son en general del tipo pistola atomizadora, las 
cuales emplean aire comprimido como fluido auxiliar. El empleo de la atomizacion 
electrostatica aumenta la confiabilidad y rentabilidad del proceso. El liquido que surge de un 
inyector se somete en este caso a la influencia de un campo electrico, obteniendose gotas 
cargadas de tamanos menores luego de la ruptura, que son guiadas hacia la superficie a 
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recubrir por fuerzas coulombianas a lo largo de las Imeas de campo electrico, reduciendo 
notablemente las perdidas de pintura. 

1.6.5 - Impresoras de chorro de tinta 

La impresion por chorro de tinta agrupa hoy en dia una gran variedad de procedimientos 
controlados electronicamente en los cuales el chorro o las gotas de tinta son depositadas sobre 
una superficie determinada. Diferentes tecnologias de impresion por chorro de tinta 
aprovechan los fenomenos electrofluidodinamicos. En esta aplicacion reviste fundamental 
importancia la supresion de gotas satelite, que alterarian la calidad del trabajo. 

1.6.6 - Produccion de polvos metalicos 

Para producir polvos metalicos puros en forma de pequenas esferas se emplea un 
metodo en el cual un cono liquido es producido por fusion en el extremo de una barra 
metalica consumible que oficia de electrodo. La aplicacion de un alto potencial de corriente 
continua en la barra permite atomizar el metal liquido en gotas esfericas que constituyen el 
polvo metalico en solidificacion. El dispositivo de atomizacion consiste en el sistema de 
alimentacion automatico de alambre metalico, la region de atomizacion y la de recoleccion 
del polvo. De esta forma pueden producirse polvos de cobre, paladio, tantalio, molibdeno y 
aleaciones, y de compuestos metalicos - ceramicos. 

1.6.7 - Atomizacion de combustibles 

El proposito de la atomizacion de combustibles en motores y quemadores es 
incrementar la relacion superficie - volumen del liquido, para de esta manera aumentar el area 
de contacto entre el combustible y el aire, favoreciendo la evaporacion del combustible y las 
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reacciones de la combustion. En forma convencional, la atomizacion se produce por medios 
mecanicos o neumaticos. EI efecto de la aplicacion de un campo electrico sobre el Hquido a la 
salida del inyector puede ser la obtencion de un aerosol de gotas mas pequenas y de tamano 
mas uniforme, lo cual favorece la combustion reduciendo las emisiones de particulas 
formadas por el combustible que no participa totalmente de la reaccion. Adicionalmente, la 
aplicacion de un campo electrico conduce al aumento del angulo del cono que forman las 
gotas, lo cual favorece el contacto entre el combustible y el aire. 

1.7 - CONCLUSIONES 

A lo largo de este capitulo se ban presentado algunas nociones elementales sobre 
chorros, tanto scan o no excitados, que serviran de base para el estudio del capitulo siguiente. 

En primer lugar se han descripto los diferentes regimenes y la variacion de la longitud 
de ruptura en funcion de la velocidad para chorros no excitados. Se ha hecho mencion a 
diferentes estudios que han establecido claramente que la ruptura de un chorro en gotas es el 
producto de la amplificacion de perturbaciones axisimetricas de longitudes de onda que estan 
dentro de un determinado rango, el cual guarda relacion con el diametro del chorro sin 
perturbar. La tasa exponencial de crecimiento de la amplitud de las perturbaciones respecto 
del tiempo, presenta un maximo en el rango antedicho. En general, las perturbaciones que 
conducen a la ruptura tienen longitudes de onda cercanas a aquella para la cual la tasa de 
crecimiento es maxima. 

Posteriormente se presentaron los chorros excitados tanto acustica como electricamente 
y los distintos procedimientos que se siguen para lograr la excitacion en cada caso. Se ha 
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hecho mencion al estudio de inestabilidades en ambos casos, y se presentaron los regimenes 
que pueden encontrarse en cada uno de ellos. 

En el caso de chorros excitados acusticamente se pusieron de manifiesto las diferencias 
que existen respecto de chorros sin excitacion en lo que a inestabilidad se refiere, y se ha 
mencionado que los modos que pueden encontrarse son los mismos que para el caso no 
excitado. 

En el caso de chorros con excitacion electrica, se han destacado las consecuencias que 
se producen en lo que respecta a la alteracion de la estabilidad del escurrimiento y la creacion 
de gotas cargadas luego de la ruptura. Esta situacion puede ser aprovechada para establecer un 
medio de control del tamano y dispersion de gotas, y del guiado de trayectoria de las mismas. 
Sin embargo el analisis de los diferentes fenomenos que ocurren en chorros electrificados no 
es sencillo debido la diversidad de regimenes, la complejidad del estudio de inestabilidades y 
la escasez de casos resueltos. Se centrara la atencion sobre ello en el capitulo siguiente, y tal 
analisis se limitara al caso particular de este trabajo, el cual presenta una relativa simplicidad. 

Finalmente se ha hecho una breve explicacion de algunas de las aplicaciones practicas 
de los chorros liquidos, las cuales unidas al interes del estudio de estos escurrimientos desde 
el punto de vista teorico (situacion esta que puede echar luz sobre otros campos de la 
Mecanica de Fluidos), hacen a la importancia del presente trabajo. 
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2.1 - EVTRODUCCION 

Habiendose presentado la dinamica de los chorros Hquidos en el capitulo anterior, en el 
presente se discuten consideraciones teoricas directamente relacionadas con el estudio 
experimental que ha sido Uevado a cabo. 

Desde un punto de vista general, en primer lugar se realiza una introduccion a la teoria 
de inestabilidades, haciendo hincapie en los diferentes tipos de analisis que pueden 
presentarse y en las conclusiones que de ellos pueden derivarse (volviendo en parte sobre 
conceptos esbozados en el capitulo anterior, e introduciendo otros). 

Posteriormente, asociado con la propagacion de pulsos en un chorro, se presentan las 
definiciones de los limites del paquete de ondas y la tasa de crecimiento, analizandose ademas 
algunas de las posibles respuestas de un chorro frente a este tipo de perturbaciones. 

Luego se centra la atencion en el estudio teorico de la influencia de los campos 
electricos sobre la propagacion de una perturbacion tipo pulso. Para ello se comienza por 
definir la geometria del problema, se establecen luego las hipotesis simplificativas, y se 
presentan las ecuaciones resultantes y las previsiones teoricas que serviran de base de 
comparacion para el estudio experimental. 

En esta circunstancia, finalmente se esta en condiciones de introducir claramente el 
objetivo del presente trabajo. 

2.2 - TEORIA DE EVESTABILmADES EN GENERAL 

En este apartado se vuelve sobre el concepto de inestabilidad presentado en el capitulo 
anterior. Se diferencian los caminos existentes para el estudio de inestabilidades, se establece 



57 



Marco teorico de las experiencias 



luego la diferencia entre las teorias lineales y las no lineales, y son introducidos los conceptos 
de inestabilidades espacial y temporal, absoluta y convectiva, y local y global. 

2.2.1 - Definicion y estudio de inestabilidades 

El estudio de la estabilidad de un sistema consiste en observar la evolucion de una 
perturbacion infinitesimal. Si esta perturbacion es amortiguada con el transcurso del tiempo el 
sistema es considerado estable (en referenda a esta perturbacion). Por el contrario si la 
perturbacion crece con el transcurso del tiempo, el sistema se considera inestable (en 
referenda a esta perturbacion) [1]. Generalizando, un sistema se considera estable si lo es en 
referenda a todas las perturbaciones posibles, y se considera inestable si lo es al menos en 
referencia a una perturbacion. 

Para efectuar el estudio de la estabilidad de un sistema existen dos caminos. El primero 
emplea los metodos globales, mientras que el otro es el analisis en modos normales. 

Los metodos globales consisten habitualmente en hallar un potencial generalizado de 
produccion de entropia y en la observacion de la estabilidad en terminos de extremos de este 
potencial. Tal potencial, si existe, es Uamado funcion de Lyapunov. Esta definicion es valida 
para metodos como el de la energia generalizada [2]. 

El analisis en modos normales consiste en descomponer la perturbacion (por ejemplo en 
serie de Fourier) y examinar la estabilidad del sistema observando cada uno de los terminos 
de la descomposicion (cada uno de los modos). 



Fisicamente, una perturbacion puede ser pensada como un "paquete de ondas" 
de extension finita, expresado por [1]: 
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(p) , > 0 ) (2-1) 

con (q) ^ < 0 ) , ^el vector posicion, t el tiempo, z la coordenada axial, 0 la 

coordenada angular, rj el valor que adopta el modo caracterizado por (k,n) y por la frecuencia 
G), para ry t dados, y tjo la amplitud infinitesimal inicial del modo normal. 

Las soluciones del tipo (2.1) podran existir solamente si <yy A: estan relacionadas 
adecuadamente. La ecuacion que determina las relaciones permisibles se conoce como 
ecuacion de dispersion [3]. En dicha ecuacion, ademas de la frecuencia coy del numero de 
onda k, se encuentran involucrados como variables algunos parametros del escurrimiento, que 
son independientes de coy k, tales como, por ejemplo, el numero de Reynolds. La ecuacion de 
dispersion se nota como sigue: 

D(co,n,k,...) = 0 (2.2) 

Las soluciones de la ecuacion de dispersion son, en el caso mas general, magnitudes 
complejas. El analisis de los valores que toman las partes real e imaginaria de las soluciones 
permite establecer en primer termino si la onda en cuestion resulta amortiguada o amplificada, 
y en segundo lugar, clasificar el tipo de inestabilidad. Sobre estos aspectos se volvera 
posteriormente a lo largo del capitulo. 

2.2.2 - Analisis lineal y no lineal 

En orden de clasificar los tipos y analisis de inestabilidades, se comenzara por 
establecer las diferencias entre las teorias lineales y aquellas no lineales. 
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Basicamente la distincion esta hecha en la manera de describir la geometria del chorro 
perturbado, cuando se descompone la perturbacion en series de Fourier. En el primer caso, de 
teorlas lineales, al desarroUar la perturbacion en series de Fourier solo se conservan los 
terminos de primer orden, mientras que en el caso de teorias no lineales, se conservan 
terminos de orden superior. 

Los analisis lineales solo son rigurosamente aplicables en tanto y en cuanto las 
amplitudes de las perturbaciones son mucho menores que el radio del chorro no perturbado, 
razon por la cual alguna de las predicciones que pueden hacerse en base a ellos, tales como la 
longitud de ruptura, solo tienen un caracter aproximado. No obstante, en un gran niimero de 
problemas estos analisis son suficientemente precisos como para dar indicaciones confiables 
sobre la estabilidad de un chorro, ya que la amplitud de las perturbaciones toma importancia 
solo en el final de la longitud intacta [4]. 

Las teorias no lineales en cambio, se ajustan mucho mejor a las observaciones 
experimentales sobre todo cerca de la zona de ruptura, en donde la amplitud de las 
perturbaciones adquiere relevancia. A pesar de que las hipotesis basicas acerca del chorro, las 
ecuaciones de movimiento y las condiciones de borde son las mismas que en las teorias 
lineales, las aproximaciones matematicas son diferentes en este caso y permiten arribar a 
conclusiones distintas [5]. Ejemplos de ello son el estudio del efecto que la amplitud inicial de 
la perturbacion tiene sobre su evolucion en el escurrimiento, y el analisis de la evolucion de 
perturbaciones que poseen fi-ecuencias armonicas de la fundamental, lo cual permite predecir 
en algunos casos, las caracteristicas de las rupturas denominadas anormales (esto es, 
formacion de gotas satelites), o bien establecer medios para suprimir las gotas satelites cuando 
las mismas no son deseadas. 
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2.2.3 - Analisis temporal y espacial 

Esta otra distincion aparece asociada al tipo de perturbacion impuesta al escurrimiento. 

En el analisis temporal se considera que la perturbacion es aplicada al flujo en una cierta 
region del espacio en el tiempo igual a cero. El desarroUo en series de Fourier de este tipo de 
perturbacion proporciona componentes con valores reales del numero de onda k, y la 
variacion en el tiempo es obtenida con las frecuencias complejas co(k) = co^(k) + icOf(k) , 

raices de la ecuacion de dispersion. En general, en este tipo de analisis, el referencial 
empleado para efectuar la descripcion matematica se mueve solidario a la perturbacion. 

En el analisis espacial se considera que la perturbacion al escurrimiento es aplicada en 
una cierta region del espacio, y que la misma sigue una cierta funcion en el tiempo. El 
desarroUo en series de Fourier de este tipo de perturbacion arroja componentes con valores 
reales de frecuencia co, y la variacion en el espacio es obtenida con los niimeros de onda 
complejos Kjp , , obtenidos de la ecuacion de dispersion. Para este caso, en 

general el referencial que se emplea es fijo al laboratorio. 

La mayoria de los trabajos de investigacion en chorros han sido dedicados a los analisis 
temporales. Ejemplos de ello y de las conclusiones a que pueden arribarse son los 
mencionados en el capitulo anterior, en ocasion de describir los chorros liquidos no excitados. 

Los analisis espaciales son generalmente empleados para describir la respuesta o 
"receptividad del flujo" para distintas frecuencias de excitacion. Un caso especial de ello lo 
constituye el denominado "problema de la senal", en donde las perturbaciones son forzadas en 
una region especifica de manera periodica, lo cual por ejemplo se corresponde con los chorros 
liquidos excitados acusticamente, descriptos con anterioridad. 
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En el caso del analisis temporal, la estabilidad del escurrimiento esta asegurada si la 
parte real de la frecuencia es negativa {co^ <0). Si la parte real de la frecuencia resulta 
positiva, el sistema es inestable {a>^ > O) . Este exponente es denominado tasa de crecimiento, 
y determina la dependencia temporal con la cual crece la amplitud del modo considerado. El 
modo que crece mas rapidamente es el mas inestable, y en general es Ilamado el modo "mas 
peligroso". La parte imaginaria de la frecuencia [co.) controla el resto de la dependencia 
temporal, y si es diferente de cero, el modo es oscilante [2]. 

Para el analisis espacial, un criterio clasico establece que un flujo es estable cuando la 
parte imaginaria del numero de onda es positiva [k^ > O) . Cuando la parte imaginaria del 

numero de onda es negativa {k. <0), el escurrimiento es inestable. Este exponente es 

denominado tasa de crecimiento e indica la dependencia espacial con la cual evoluciona la 
perturbacion. 

En el analisis espacial, aunque las componentes de la perturbacion no posean 
frecuencias imaginarias, un segundo criterio mas general que permite determinar si la onda es 
amplificada (lo que se corresponde con el caso anterior de inestabilidad) o no transmitida 
(situacion que se corresponde con el caso de anterior de estabilidad), requiere considerar 
como complejos a la frecuencia y al numero de onda [6]. Una onda de valor complejo k{(Dg) 

para una frecuencia real dada (Og , sera amplificada si la funcion k^((o) cambia de signo segun 

0)^ vaya de +oo a 0 (con m = m^ +i(Oi, y <y^ =g>q). Si ki((o) no cambia de signo la onda 

entonces no se transmite. Por supuesto, este criterio proporciona la misma informacion que el 
criterio anterior. 



62 



Marco teorico de las experiencias 



En general, los analisis temporal y espacial arrojan la misma informacion solo cuando 
se cumplen las condiciones especificadas en el teorema de Gaster [7]. 

2.2.4 - Inestabilidades absoluta y convectiva 

No hay una regla absoluta para determinar de manera anticipada cual analisis es el 
apropiado para un problema en particular. La razon es que las teorias espaciales describen la 
respuesta de un sistema forzado solamente en el limite del tiempo tendiendo a infinito. Para 
decidir en que manera un sistema evoluciona (temporalmente, espacialmente, u otro en el cual 
tanto k como co sean complejos), es necesario evaluar la respuesta asintotica del mismo ante 
una perturbacion tipo pulso, considerando como complejos tanto al numero de onda como a la 
frecuencia [8]. Esta respuesta esta asociada con la presente clasificacion, segiin se vera. 

Toda perturbacion puede ser pensada como un paquete de ondas en el espacio, de 
dimension finita, como se ha destacado con anterioridad. Para un escurrimiento inestable, 
pueden distinguirse dos tipos de inestabilidad, la convectiva y la absoluta [6]. 

Si el paquete de ondas creciente es arrastrado alejandose de la fuente de perturbacion, la 
inestabilidad es convectiva. Esto finalmente produce la ruptura del chorro cuando la amplitud 
es lo suficientemente grande, hecho que toma lugar en una posicion alejada de la fuente de 
perturbacion (en el contexto matematico puede decirse que si no se Uegara a la ruptura, en 
todo punto de la superficie del chorro la amplitud de la perturbacion tiende a cero cuando el 
tiempo tiende a infinito). 

Por el contrario, en un flujo inestable absolutamente, una perturbacion localizada se 
extiende aguas arriba y abajo de la fiiente de perturbacion, contaminando el escurrimiento por 
completo y creciendo indefinidamente en todo punto del espacio (en terminos matematicos. 
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puede decirse que si no se Uegara a la ruptura, en todo punto del flujo la amplitud de la 

perturbacion tenderia a infinito cuando el tiempo tiende a infinito). 

El paquete de ondas que representa a la perturbacion, se propaga en el escurrimiento 
con una cierta velocidad de grupo, que esta definida mediante: 

' dD/dco dk dk 

Algunos autores asocian el caracter absolute o convectivo de la inestabilidad a flujos 
subcriticos o supercriticos, en los cuales la velocidad del escurrimiento es menor o mayor que 
la velocidad de grupo, respectivamente. De acuerdo con [6], una inestabilidad absoluta para 
un cierto referencial dado (por ejemplo uno que se desplaza a la velocidad del chorro), puede 
transformarse en una inestabilidad convectiva en otro referencial (por ejemplo uno fijo en el 
origen del chorro). 

Asociado con ello aparece un primer criterio para diferenciar una inestabilidad absoluta 
de una convectiva. El mismo establece que, para un referencial dado, la inestabilidad es 
absoluta cuando la velocidad de grupo de la perturbacion es nula. 

Un segundo criterio [6, 9], se basa en considerar como complejos tanto al numero de 
onda como a la frecuencia. En este caso puede fijarse un valor para co^ y determinar los 
valores de k, para un cierto valor de a)^ , que puede ser variado. La ecuacion que permite la 
determinacion de k arroja una raiz doble que esta representada por dos puntos en el piano 
complejo k. Los diferentes valores de oj^cambian la posicion de los puntos k solucion de la 

ecuacion. Puede existir entonces una cierta combinacion de co^y co^, identificada como 
Q}^ =(jq,^+cq^ (punto comunmente denominado frecuencia absoluta), para la cual la raiz k{cL)) es 
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unica, es decir que las dos raices se confunden. El criterio establece que un sistema es 
absolutamente inestable si co^^ > 0 , y convectivamente inestable si co^^ < 0 (la parte 
imaginaria de la frecuencia absoluta se denomina tasa de crecimiento absoluta). Cuando 
(o.^ — > +00 , las dos raices se situan de uno y otro lado del eje real del piano complejo k. 

Cuando se da la situacion de igualdad de las dos raices, se cumple la condicion de nulidad de 
la velocidad de grupo, en concordancia con el criterio anteriormente expuesto. 

Para el segundo de estos criterios, Gaster [10] sugiere considerar la respuesta forzada 
del sistema ante una perturbacion tipo pulso, ya que tal perturbacion contiene modos 
caracterizados por todas las frecuencias y niimeros de onda, y de esa manera cualquier 
singularidad se vera reflejada en la respuesta del sistema. Generalmente, en este tipo de 
analisis solo se toma en consideracion a la singularidad que describe el modo de mayor tasa 
de amplificacion. 

La mayor parte de la literatura ha sido consagrada al estudio de inestabilidades 
convectivas. En el campo de las inestabilidades absolutas, Leib y Goldstein [11,12] han 
establecido que si se desprecian efectos gravitatorios y de la atmosfera circundante, un chorro 
liquido con perfil inicial de velocidad parabolico, el cual se relaja hacia otro tipo tapon, y con 
un niimero de Reynolds bajo, puede volverse absolutamente inestable por accion de las 
fuerzas de tension superficial. En sus trabajos se muestran curvas en las cuales se representa 
el pasaje de uno a otro tipo de inestabilidad, en funcion de un parametro caracteristico del 
perfil inicial de velocidad y del numero adimensional de Weber, el cual representa la relacion 
entre las fuerzas de inercia y las de tension superficial, y esta definido por: 

We=''^^ (2.4) 

r 
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con p la densidad del Hquido, Uo la velocidad media del escurrimiento, a el radio del 
chorro y / la tension superficial. Para el caso de un chorro con perfil de velocidad inicial tipo 
tapon, muestran que el valor de We para el cual la inestabilidad cambia entre absoluta y 
convectiva, tiene un valor asintotico de aproximadamente 3,1 segun Re tiende a infinito. 
Respecto del efecto de la viscosidad del liquido, concluyeron que cuando la misma aumenta y 
para Re bajos, disminuye el We critico que marca el cambio entre inestabilidad absoluta y 
convectiva. En cualquiera de los dos casos, el aumento de We para valores determinados de 
los parametros, hace que se vea reducida la tasa de crecimiento absoluta, o en otra palabras, 
que la velocidad de propagacion de las ondas aumente y el flujo "tienda a ser mas 
convectivo". Con este tipo de analisis, puede intentarse una explicacion de la transicion entre 
el regimen de goteo y el de un chorro completamente establecido, cuando se incrementa la 
velocidad del escurrimiento. 

Lin y Lian determinaron que la inestabilidad absoluta en chorros liquidos tambien 
puede darse cuando se considera el efecto de un medio gaseoso circundante, tanto sea para un 
gas viscoso [13], como no viscoso [14]. Para el caso del gas no viscoso, encontraron que el 
aumento de la relacion entre las densidades del gas y del liquido, ateniia el efecto antes 
mencionado de la viscosidad del liquido. Con el aumento de la relacion de densidades, la tasa 
de crecimiento absoluta tiende a incrementarse, la velocidad de propagacion de las ondas se 
ve disminuida, por lo que el escurrimiento "tiende a ser mas absoluto". Cuando el efecto 
viscoso del gas es considerado, concluyeron que el dominio de inestabilidad absoluta resulta 
ampliado significativamente. 

Zhou y Lin [15], encontraron que el fenomeno de inestabilidad absoluta en chorros de 
fluidos no viscosos y compresibles, que se desarroUan dentro de otro fluido no viscoso y 
compresible, aparece no solo asociada a efectos de la tension superficial, sino tambien a 
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efectos de fluctuacion en la presion interfacial. La inestabilidad absoluta resulta de esta 
manera promovida por la compresibilidad del fluido circundante. El caracter absolute de la 
inestabilidad puede verse reducido, en el caso analizado, tanto sea mediante la reduccion de la 
tension superficial como por el incremento de la compresibilidad del fluido del chorro. 

En el campo de chorros liquidos excitados electricamente, Artana et al [16] 
concluyeron que el efecto de un campo electrico coaxial con el escurrimiento, para un chorro 
liquido de alta conductividad electrica, es similar al mostrado en [14] cuando aumenta la 
relacion de densidades. Un aumento en la intensidad del campo electrico puede cambiar el 
caracter de la inestabilidad, de convectiva a absoluta, disminuyendo la velocidad de 
propagacion del paquete de ondas, pudiendo el flujo pasar de supercritico a subcritico. 

Un estudio experimental sobre inestabilidades absolutas en chorros liquidos fue 
efectuado por Vihinen et al [17]. EUos muestran en su trabajo la transicion de uno a otro tipo 
de inestabilidad en funcion de parametros adimensionales del escurrimiento, los numeros de 
Reynolds y Weber. La conclusion a la que arribaron es que las predicciones teoricas 
concuerdan cualitativamente con sus resultados experimentales, en el rango de Re y We 
estudiado y pasible de ser comparado con resultados teoricos (para escurrimientos con valores 
muy bajos de Re y We, las teorias no son precisas debido a dificultades de indole 
computacional). 

Un chorro convectivamente inestable de aceite SAE 10 es mostrado en la figura 2.1. 
Puede apreciarse que se trata de un chorro establecido, perturbado axisimetricamente (modo 
n=0), en el cual las perturbaciones se alejan de la fuente y fmalmente provocan la ruptura en 
gotas. 
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Figura 2.1 - Chorro convectivamente inestable. Re=5,4; We=l 7,857 [17] . 
La figura 2.2 es un ejemplo de inestabilidad absoluta, en un chorro de glicerina. Luego 
de la aparicion de la inestabilidad absoluta, la perturbacion se propaga aguas arriba de forma 
repentina agolpandose hacia el orificio y fDrmando un pendiente, mientras que se propaga 
aguas abajo creciendo lentamente sobre el fino chorro que emana del pendiente. Entre 0,2 y 
0,4 segundos luego de que la gota se ha separado, el filamento liquido resulta estrechado en el 
pendiente, debido a la accion de la tension superficial, hasta que se produce su separacion. 
Tambien puede observarse que la amplitud de la perturbacion crece en el tiempo, y que el 
pendiente va tomando una forma cada vez mas esferica, en tanto que el flujo por el orificio 
permanece constante. El modo de inestabilidad observado en esta figura se corresponde con el 
modo global que se presentara posteriormente (luego de la aparicion de la inestabilidad 
absoluta, la amplitud de las perturbaciones crece a lo largo de toda la superficie del chorro). 

Observando ambas figuras, puede verse que en el caso del problema de la serial, un 
sistema inestable solamente puede comportarse como convectivo. Esto se debe a que en el 
caso de inestabilidad absoluta la perturbacion crece indefinidamente con el tiempo en todo 
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punto del escurrimiento, de manera tal que el establecimiento de un patron de ondas en 
regimen permanente es imposible [6]. 




Figura 2.2 - Chorro absolutamente inestable. Re=0,082; We=0,349. Izquierda: 0,2 seg. 
despues de la separacion de la gota. Derecha: 0,4 seg. despues del goteo [17]. 

En resumen, la inestabilidad de un escurrimiento puede ser clasificada en convectiva o 

absoluta de acuerdo a si el paquete de ondas representativo de las perturbaciones se propague 

aguas abajo de la fuente o en ambas direcciones, respectivamente. El pasaje de uno a otro tipo 

de inestabilidad ha sido estudiado teoricamente, encontrandose las relaciones apropiadas entre 

los parametros del escurrimiento para que este hecho tome lugar. Los efectos producidos por 

la viscosidad del chorro, la atmosfera circundante y la excitacion electrica entre otros, sobre 

los parametros antes mencionados, tambien han sido investigados en forma teorica. En 

cambio en el campo experimental, la bibliografia sobre este fenomeno es mucho mas escasa 

hasta el momento. 
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2.2.5 - Inestabilidades local y global 

Los analisis temporales y espaciales discutidos anteriormente, por lo general estan 
restringidos a la estabilidad local del perfil de velocidades en un sitio tipico del escurrimiento, 
el cual es invariante en la direccion de la propagacion de las perturbaciones. 

Sin embargo, cuando el escurrimiento basico difiere de un sitio a otro de manera 
significativa, el analisis de la estabilidad requiere considerar no solamente la estabilidad del 
perfil local de velocidad en un sitio caracteristico (lo cual constituye un analisis local), sino la 
estabilidad del escurrimiento por completo (que es el caso de un analisis global). 

El analisis de la estabilidad de chorros en donde se toman en cuenta efectos 
gravitatorios, la relajacion del perfil inicial de velocidad, o la accion viscosa de la atmosfera 
circundante, entre otros fenomenos, requiere de los metodos globales ya que en esos casos se 
presenta variacion del flujo basico de un lugar a otro. 

Algunos autores [8, 18, 19], sugieren que en flujos de superficies libres que se 
desarroUan espacialmente, existe una vinculacion entre las propiedades de inestabilidades 
locales y globales, si en la escala de la longitud de onda de la perturbacion, el flujo basico se 
desarroUa lentamente en la direccion del escurrimiento. Tal relacion puede ser resumida como 
sigue: la inestabilidad global puede ser promovida por la existencia de inestabilidad local y 
absoluta, en tanto que la ausencia de inestabilidad local y absoluta es condicion suficiente 
para asegurar la estabilidad global cuando los efectos de borde y de realimentacion son 
despreciables. Esta vinculacion provee un medio de diagnostico de relativa sencillez para 
flujos globalmente inestables. 

Analogamente a la clasiflcacion de inestabilidades absoluta o convectiva, para este tipo 
de escurrimiento establecieron una distincion para el comportamiento global entre 
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"amplificador" y "oscilador". El primer caso corresponde a la situacion en la cual el sistema 
retoma al estado no perturbado cuando cesa la excitacion externa, comportandose como un 
amplificador que selectivamente amplifica ruidos extemos. El segundo caso se refiere a la 
situacion en donde el sistema, partiendo del estado no perturbado, desarrolla oscilaciones 
autoexcitadas que crecen con el tiempo en cualquier lugar fijo, al ser perturbado mediante un 
pequeno impulso. En este ultimo caso, el sistema muestra una respuesta global intrinseca que 
se desarrolla con frecuencias caracteristicas propias del sistema oscilatorio. En general, tal 
respuesta es independiente de la naturaleza del impulso inicial, y eventualmente de cualquier 
accion adicional externa, forzada y continua, de bajo nivel. 

Estos autores han comprobado la validez de sus analisis mediante estudios 
experimentales que incluyen, entre otros, estelas de vortices detras de cuerpos romos 
bidimensionales, transicion al regimen de goteo en chorros liquidos y chorros calefaccionados 
de baja densidad. 

Un estudio experimental sobre chorros en los cuales el perfil de velocidad se relaja 
desde uno parabolico en el inicio hacia otro tipo tapon fue efectuado por Debler y Yu [20]. 
EUos concluyeron que la teoria desarroUada en [11], coincide cualitativamente con sus 
resultados en relacion a la tasa de crecimiento de las perturbaciones, mostrando que la 
relajacion del perfil de velocidad conduce a una reduccion del caracter inestable del flujo. En 
tanto que ellos determinaron una relacion de 2:1 entre las tasas de crecimiento de un chorro 
con perfil de velocidades piano y uno con perfil parabolico en el inicio, en [11] se muestra 
una relacion de 4:1. En su trabajo ademas proponen un modelo para este tipo de problemas, el 
cual muestra concordancia con las observaciones experimentales. 

Yakubenko [21] efectuo un analisis lineal de la estabilidad global en un chorro liquido 
que emana de un inyector inclinado respecto al campo gravitatorio. La conclusion a la que 
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arribo en su trabajo es que un modo global autoexcitado aparece aun cuando en todo punto del 
escurrimiento, el flujo es local y convectivamente inestable, esto es cuando no hay regiones 
de inestabilidad local y absoluta, y que el inyector juega un rol preponderante en la aparicion 
del modo global al influir en el perfil inicial de velocidad en el chorro. Sus conclusiones 
teoricas ademas aparecen apoyadas por un sencillo estudio experimental incluido en el mismo 
trabajo. El modo global esta compuesto de tres ondas, las cuales interactiian en varios puntos 
del escurrimiento, cercanos al orificio. El resultado de este trabajo es sorprendente, en el 
sentido de que contradice la conexion entre inestabilidades locales y globales comentada 
anteriormente. 

Mediante un analisis lineal, Yakubenko [22] encontro que un chorro circular no 
viscoso de gran longitud, pero finita, se comporta inestable globalmente cuando son fijadas 
las condiciones de borde adecuadas en el comienzo y en el final del escurrimiento. La 
situacion de inestabilidad global se presenta aun cuando el correspondiente flujo infinito es 
convectivamente inestable. El modo global esta compuesto por un par de ondas que se 
propagan en direcciones opuestas y que se reflejan una hacia otra en los hordes del 
escurrimiento. 

Para un chorro en caida libre, Le Dizes [23] desarroUo un modelo lineal de estabilidad 
global, asumiendo que las fuerzas de inercia y las de tension superficial gobiernan la dinamica 
del chorro por sobre las gravitacionales, despreciando el efecto de la densidad del gas 
circundante y asumiendo que el perfil de velocidades es piano en cada seccion. Con tal 
analisis encontro que un modo globalmente inestable, compuesto de tres ondas locales planas 
autoexcitadas, y que interactuan en el orificio donde es formado el chorro, se presenta por 
encima de la transicion entre las inestabilidades absoluta y convectiva locales. Estos 
resultados concuerdan cualitativamente con los mostrados en [8, 18]. Cuando el modo global 
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aparece, el flujo se comporta como un oscilador, en donde fundamentalmente las 
perturbaciones crecen modificando el flujo basico cerca del orificio. Esta modificacion es 
importante para explicar la transicion al regimen de goteo, cuando no intervienen efectos no 
lineales que tienden a estabilizar los modos globales. 

Artana y Seguin [24] analizaron el caso de un chorro en caida libre y excitado 
electrostaticamente. Suponiendo a la evolucion del chorro gobemada fundamentalmente por 
la accion de la gravedad, establecieron las relaciones entre las fuerzas de tension superficial, 
de inercia, y electricas, que aseguran la supresion de modos globales en chorros no viscosos y 
de alta conductividad electrica. En su trabajo muestran que la presencia de inestabilidad local 
absoluta en el escurrimiento, promueve la aparicion de modos globales en las cercanias del 
orificio en donde se forma el chorro, cuando existe un fiaerte campo electrico. En general, 
sugieren que chorros de baja velocidad excitados con campos electricos intensos, son 
escurrimientos en donde la estabilidad global debe ser cuidadosamente estudiada. 

Resumiendo, otro tipo de inestabilidad, la global, puede aparecer cuando por distintas 
razones, el flujo basico muestra diferencias en el comportamiento de un lugar a otro del 
escurrimiento. Cuando un flujo es globaknente inestable, puede comportarse como un 
amplificador de ruidos extemos, o como un oscilador que muestra ondas crecientes con el 
tiempo en un lugar fijo, que son autoexcitadas, y que poseen fi-ecuencias caracteristicas. La 
conexion entre las propiedades de inestabilidad local y el comportamiento global de un flujo 
no esta clara y es un motivo de controversias. No obstante, algunas observaciones sostienen la 
idea de que la existencia de inestabilidades locales y absolutas en un escurrimiento, es una 
condicion necesaria para la aparicion de modos globales. La presentacion de algunos trabajos 
en el campo de inestabilidades globales, hecha precedentemente, es una muestra de ello. 
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2.3 - PROPAGACION DE UN PULSO 

La respuesta de un flujo ante una perturbacion tipo pulso impuesta al mismo, ha sido 
estudiada por Gaster [10], entre otros [8, 18]. El problema considerado es aquel de un flujo 
paralelo (pudiendo aplicarse sus conclusiones a un chorro liquido circular), excitado con una 
perturbacion pulsante en el espacio y en el tiempo. Ante tal perturbacion, el sistema responde 
desarrollando una onda viajera bidimensional, la cual para el case de un chorro circular de 
radio r y coordenada axial z, muestra una dependencia espacio-temporal que puede ser 
descripta mediante la siguiente ecuacion, con ky co magnitudes complejas, siendo t el tiempo 
transcurrido luego de la aplicacion de la perturbacion: 



De particular interes resulta el comportamiento asintotico de la respuesta para ? — > oo , 
ya que esto provee un medio de analisis del escurrimiento en terminos de inestabilidades 
absoluta o convectiva. La respuesta asintotica, o funcion de Green, es obtenida 
convenientemente en terminos de A: y <y, en el espacio de Fourier. Esta funcion esta definida 
como: 



Evaluando la transformada inversa de Fourier, se halla que la respuesta sobre cada rayo 
z/t=cte., esta gobemada por el modo menos estable, con una velocidad de grupo real igual a 

z/t. Para este modo, los correspondientes niimero de onda k y frecuencia co , relacionados a 



y/{r,z,t) = (l){r)-e'^''-'-'-''^ 



(2.5) 




(2.6) 
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traves de la ecuacion de dispersion, satisfacen entonces las siguientes ecuaciones (que surgen 
a partir de (2.3): 



dco^ 



(r)=^ (2.7) 



Sobre cada rayo z/t, el modo experimenta un crecimiento en el espacio y en el tiempo, 
con una tasa definida como sigue (en donde el miembro de la derecha esta escrito en la fonna 
de un crecimiento temporal): 



CT^=-k;z + 0),(k^)t = - 



k*i ■--C0j[k*) 



(2.9) 



Si un escurrimiento es inestable, el paquete de ondas resulta confinado entre dos rayos 
para los cuales <T^ = 0 (los rayos de "cabeza" y de "cola"), y la maxima tasa de crecimiento 
asociada con el modo mas amplificado, aparece para un rayo z/t determinado, para el cual 



(7=0) 

I I. max 



La figura 2.3 es un ejemplo de posibles respuestas ante una perturbacion pulsante, para 
diferentes rangos de un parametro de control R generico. Las perturbaciones crecen 
exponencialmente sobre todos los rayos contenidos dentro de los limites indicados, y son 
amortiguadas fuera de estos limites. Un primer valor critico Rc marca la transicion entre un 
flujo estable representado por (a), a uno convectivamente inestable mostrado en (b). Existe 
otro punto de transicion dado por Rt > Rc, por encima del cual el flujo se vuelve 
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absolutamente inestable (c). Puede observarse en este caso que la perturbacion se propaga 
tambien aguas arriba de la fuente. 



I 



/ 




8 <Rc 

Resp. Estable Resp. Convectiva Resp. Absoluta 

Figura 2.3 - Representaciones de respuestas impulso Hpicas [18]. 
Estas representaciones se refieren a un flujo que se comporta de manera localmente 
inestable. Las distintas respuestas para flujos que muestran modos de inestabilidad global, no 
seran analizadas aqui pero pueden ser consultadas en [18]. 

2.4 - EFECTO DE UN CAMPO ELECTRICO EN LA RESPUESTA DE UN CHORRO 
LIQUIDO FRENTE A UNA PERTURBACION TIPO PULSO 



Se ha mencionado que en un escurrimiento inestable, la respuesta a una perturbacion 
pulsante puede ser caracterizada por dos parametros: la tasa de crecimiento de la perturbacion 
y la velocidad con la cual la misma se propaga en el flujo. Respecto del ultimo parametro, se 
vio que un campo electrico reduce la velocidad de propagacion haciendo en teoria posible 
gobemar el caracter absoluto o convectivo de la inestabilidad de un chorro [16]. No obstante 
ello, con tal analisis no pueden especificarse para cada posicion axial del escurrimiento tanto 
la tasa de crecimiento como la expansion del paquete de ondas. 

Este apartado esta consagrado entonces a establecer la influencia de un campo electrico 
sobre los limites del paquete de ondas y la tasa de crecimiento de perturbaciones pulsantes 
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para diferentes posiciones axiales, en funcion del tiempo, siguiendo el trabajo de Artana y 
Seguin [25]. Los resultados surgen de un estudio de inestabilidades lineal, local y en modos 
normales. 



2.4.1 - Geometria del problema e hipotesis simplificativas 

En la figura a continuacion se representa la situacion no perturbada del problema en 
estudio: un cilindro liquido infinito de radio a escurriendo a velocidad uniforme Og, en un 

medio circundante de aire sujeto a presion atmosferica, en el interior de un electrodo 
cilindrico coaxial de radio b sometido a un potencial electrico diferente al del liquido. 
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Figura 2.4- Esquema del problema analizado. 
Se consideran las siguientes hipotesis simplificativas: 

• El chorro se supone no viscoso, incompresible e isotermico. 

• Se admite que no hay transferencia de masa entre el chorro y el medio circundante. 

• Se desprecian los efectos de los campos gravitacional y magnetico. 

• El perfil de velocidad constante es valido para todo el escurrimiento, es decir no 
tiene dependencia de la coordenada axial. 



77 



Marco teorico de las experiencias 



• Se acepta que el tiempo de relajacion de cargas del liquido es muy corto, entonces la 
superficie del chorro es equipotencial. 

2.4.2 - Ecuacion de dispersion y sistema de ecuaciones del paquete de ondas 

Para la situacion no perturbada y en base a las hipotesis anteriores, los campos de 
velocidad y electrico se describen mediante las siguientes expresiones, en donde el subindice 
1 indica la fase liquida y el 2 la fase gaseosa: 



con (jg la densidad superficial de carga y Sg la constante dielectrica del vacio. 

Considerando una perturbacion tipo pulso en la superficie del chorro, denotada por 
=a + S(t)S(r) , la respuesta del sistema ante la misma ha sido establecida en 2.3. En el 

caso de un chorro electrificado, la diferencia con lo expuesto en aquella ocasion aparece en la 
expresion de la ecuacion de dispersion, la cual en este caso toma la forma: 



en la cual coy k son la frecuencia y el numero de onda adimensionales, Eu es el numero 
adimensional de Euler, (que relaciona las fuerzas de tension superficial con las fuerzas de 
inercia) y Eug es el numero adimensional de Euler electrico (que relaciona las fuerzas 
electricas con las fuerzas de inercia). Estos adimensionales estan definidos mediante: 



£,=0 



a,=o 



2n 



CTg-a 



D( CO, k, Eu, Eue ) = 0 



(2.10) 



co = co 



Uo 



a 



k = ak 



Eu = 



r _ 1 

PjUga We 



Eu^ = 
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donde es la densidad del liquido y /es la tension superficial. 

Detalles de la obtencion de la ecuacion de dispersion y su expresion completa pueden 
verse en [26, 27]. Respecto de su deduccion, aqui solo se puntualizara que como el paquete de 
ondas modifica la geometria del sistema, y provoca las perturbaciones sobre los campos de 
velocidades y electrico, las ecuaciones de balance del escurrimiento perturbado no seran las 
mismas que las correspondientes a la situacion no perturbada. Para la deduccion de la 
ecuacion de dispersion se requiere entonces hallar los campos de velocidades y electrico en la 
situacion perturbada y el posterior planteo de la conservacion de la cantidad de movimiento 
en la interfaz. 

La ecuacion empleada en este estudio surge de despreciar el efecto de la atmosfera 
circundante y de considerar una perturbacion axisimetrica (esto es considerar solamente el 

modo n=0, y la relacion de densidades = — = 0), ademas de suponer que el radio del 

Pi 

electrode es mucho mayor que el correspondiente al chorro {b»a). Con tales 
simplificaciones, la ecuacion de dispersion resulta: 

Eu-k-(l-k')+A-{k-coy -Eu^ ■k-{l + k-g)=0 (2.11) 

con 

A = W a = ^^ c = - 
I,(k) ' K„(k) ' a 

en donde Ig.Kg son las funciones de Bessel modificadas de primera y segunda especie, 
y ' indica la derivada rerspecto del argumento. 
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Para cada rayo z/t, la respuesta del sistema ante la perturbacion esta dominada por el 
modo menos estable, caracterizado por k y . La resolucion de este problema requiere 
considerar las condiciones sobre el grupo de velocidad presentadas en las ecuaciones (2. 7) y 
(2.8), y la separacion de la ecuacion de dispersion (2.11) en sus partes real e imaginaria, lo 
cual conduce al siguiente sistema de ecuaciones: 



„ dD/dk da)/k*) z 

Re = - = — 

dD/dco dk t 

j^^Diei^Jjhm^O (2.12) 
dD/dco dk 

ReD = 0 
ImD = 0 



Representando entonces la tasa de crecimiento en funcion de la relacion z/t es posible 
obtener los limites del paquete de ondas, que estan definidos por los rayos z/t para los cuales 
la amplificacion es nula. El paragrafo siguiente esta consagrado a tal analisis. 

2.4.3 - Previsiones teoricas 

El sistema de ecuaciones (2.12) ha sido resuelto numericamente para distintos 
escurrimientos inestables. En este apartado se muestran los resultados de tal analisis en forma 
grafica y se discuten los efectos de los parametros del flujo sobre la inestabilidad del sistema. 

En la figura 2.5 se muestra la tasa de crecimiento (adimensionalizada con el radio y la 
velocidad del chorro) en funcion de z/t (adimensionalizado con la velocidad del chorro), para 
diferentes valores del numero de Weber. Puede notarse que el incremento del valor numero de 
Weber, hace que los limites del paquete de ondas tienden a acercarse, resultando de tal forma 
mas compacto, y que la maxima tasa de crecimiento corresponde al rayo z/t=l. 
Adicionalmente, puede observarse una simetria respecto del rayo z/t=l, de los valores de z/t 
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que anulan la tasa de crecimiento. La propagacion del paquete de ondas aguas arriba de la 
fuente pulsante, ocurre cuando el rayo z/t=0 tiene una tasa de crecimiento positiva, lo cual 
corresponde a un flujo absolutamente inestable, y esto toma lugar para bajas velocidades 
correspondientes a We=3,l en concordancia con [11]. Analizando esta figura puede 
concluirse que las fuerzas de inercia tienden a estabilizar el escurrimiento ante una 
perturbacion pulsante, y promueven la propagacion de un paquete de ondas mas compacto. En 
cambio, para la tension superficial, el efecto es opuesto al anterior. 
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Figura 2.5 - Tasa de crecimiento adimensional en funcion de z/ty con We como pardmetro. 

La figura 2.6 muestra la tasa de crecimiento adimensional en funcion de la relacion 
adimensional z/t, para diferentes valores del numero de Euler electrico y para un valor dado 
de We. Puede observarse que la aplicacion de un campo electrico, no tiene efecto sobre la 
asimetria del pulso respecto del rayo adimensional z/t=l anteriormente citada, pero el paquete 
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de ondas puede propagarse aguas arriba si el campo electrico es lo suflcientemente intense. 
Valores elevados de Eue (campos electricos intensos), promueven la propagacion de un pulso 
mas expandido, ya que los limites del paquete de ondas resultan mas separados. En 
conclusion, el efecto de un campo electrico es opuesto al de las fuerzas de inercia ya que tiene 
un efecto desestabilizante sobre una perturbacion tipo pulso, incrementando la tasa de 
crecimiento segun la cual la perturbacion se amplifica en cada rayo z/t, en tanto que expande 
el paquete de ondas. Esta accion del campo electrico se ve reducida para altos valores del 
niimero de Weber (o para chorros de alta velocidad). 
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Figura 2.6- Tasa de crecimiento adimensional en funcion de z/ty con Eue como pardmetro, 

para We=5. 

En la figura 2.7 se muestra la tasa de crecimiento adimensional correspondiente al rayo 
adimensional z/t=l, en flincion de Eue, en las mismas condiciones de flujo que en la figura 
anterior. 
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Figura 2.7 — Maxima tasa de crecimiento adimensional en funcion de Rae, para We=5. 



A titulo de ejemplo de las previsiones teoricas acerca de los limites del paquete de 
ondas, se presenta en la figura 2.8 un grafico de los valores que adoptan las magnitudes 
adimensionales z/ti y z/tn, que identifican a los rayos de cabeza y de cola del paquete de 
ondas, en funcion de Eue para We=5. En adelanto a la presentacion del estudio experimental, 
puede decirse aqui que para un escurrimiento de velocidad Uo y una posicion z dadas, la 
extension del paquete de ondas en el tiempo puede calcularse con la expresion (2.14). Esto 
servira para la comparacion con los tiempos medidos sobre los pulsos. 



At = lzlt)r-{zlt\-'\^ (2.14) 
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Figura 2.8- z/ti y z/tn en funcion de Eue para We=5. 



2.5 - OBJETIVO DEL TRABAJO 



0,75 



El objetivo de este trabajo, es analizar experimentalmente la influencia de un campo 
electrico, sobre los Hmites y la tasa temporal de crecimiento del paquete de ondas resultante 
de una perturbacion tipo pulso sobre la superficie de un chorro Kquido de baja velocidad, 
convectivamente inestable, para diferentes posiciones segun la coordenada axial del chorro. 

Los resultados de este estudio experimental, proporcionaran el marco comparativo 
necesario para establecer conclusiones concernientes a la validez de la teoria expuesta en 2.4. 
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2.6 - CONCLUSIONES 

Este capitulo ha sido consagrado a la presentacion del estudio teorico de la influencia de 
un campo electrico sobre la propagacion de una perturbacion tipo pulso, en un chorro Hquido. 

Con anterioridad a ello, se han discutido otros aspectos teoricos de la inestabilidad de 
chorros, con la intencion de lograr una cabal comprension del fenomeno que aqui se trata. 

El concepto de inestabilidad que habia sido desarroUado en el capitulo anterior ha sido 
revisado, incluyendose la diferenciacion entre las vias usuales de enfrentar el modelado 
matematico de este tipo de problemas. Para el caso del metodo en modos normales, la 
importancia de la ecuacion de dispersion ha sido destacada, en tanto que fija las relaciones 
permisibles entre la frecuencia y el numero de onda, dando lugar a diferenciar entre 
perturbaciones que se amplifican o amortiguan, y a la clasificacion de la inestabilidad. 

Con posterioridad se trataron las diferencias existentes entre las teorias lineales y no 
lineales, y entre las inestabilidades temporales y espaciales, volviendo sobre algunas 
consideraciones presentadas en el capitulo anterior, al tratar los chorros sin excitacion y con 
excitacion acustica. 

Los conceptos de inestabilidades absoluta y convectiva, y local y global, fueron 
presentados a continuacion. La influencia que tienen los parametros de los escurrimientos 
sobre la aparicion y evolucion de estas inestabilidades, como asi tambien las conclusiones a 
las que pueden arribarse con estos tipos de analisis, han sido destacadas en ocasion de 
presentar en forma resumida algunos trabajos de investigacion sobre estos campos. 

Luego, la respuesta asintotica de un chorro frente a una perturbacion tipo pulso ha sido 
analizada. Se ha hecho hincapie sobre dos aspectos de fundamental importancia para el 
presente trabajo, como son las definiciones de los limites y de la tasa de crecimiento del 



85 



Marco teorico de las experiencias 



paquete de ondas resultantes de la perturbacion. Apoyando a la definicion de la frontera del 
paquete de ondas, ha sido incluida una representacion grafica de algunas de las posibles 
respuestas de un chorro ante este tipo de perturbacion, que tienen implicancias sobre este 
trabajo. 

Finalmente, se Uega al desarroUo del estudio teorico de la manera en la cual un campo 
electrico, actuando sobre la superficie del chorro, afecta a la respuesta aludida anteriormente. 
La geometria del problema, las hipotesis simplificativas, y las ecuaciones resultantes del 
estudio de la inestabilidad del escurrimiento mediante un analisis lineal, local y en modos 
normales, han sido presentadas. 

De la observacion de los resultados obtenidos, surge que el efecto de un campo electrico 
sobre la propagacion de una perturbacion tipo pulso en un chorro liquido, es diferente a la 
accion de las fuerzas de inercia y similar al efecto de las fuerzas de tension superficial. Un 
incremento de las fuerzas de inercia tiende a estabilizar una perturbacion pulsante y promueve 
la propagacion de un paquete de ondas mas compacto, mientras que el incremento de las 
fuerzas de tension superficial tiene un efecto contrario al anterior. Los resultados expuestos 
muestran que el efecto de un campo electrico actuando en la superficie del chorro es 
incrementar la separacion entre las fronteras del paquete de ondas (la cual puede darse aiin 
aguas arribas de la fiiente de perturbacion para campos suficientemente intensos), y ademas 
incrementar la tasa de crecimiento con la cual la perturbacion es amplificada en cada rayo 
z/t=cte. Esta accion se ve reducida cuando el numero adimensional de Weber adopta un valor 
alto (o para velocidades altas del chorro). 
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El estudio experimental de este fenomeno es el objeto del presente trabajo. Para dar 
cumplimiento a ello un dispositive experimental debe ser desarrollado y la presentacion de 
sus caracteristicas sera el tema a tratar en el siguiente capitulo. 
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3 - MONTAJE EXPERIMENTAL 
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3.1 - EVTRODUCCION 

En este capitulo se describe el montaje experimental disenado y construido para realizar 
las mediciones. 

Las experiencias pueden ser realizadas siguiendo alguno de los dos metodos siguientes: 
analizando las perturbaciones en el chorro mediante imagenes obtenidas en diferentes 
regiones del mismo, quedando fija la variable tiempo, o bien estudiando la evolucion del radio 
del chorro como funcion del tiempo, para distintas posiciones fijas segiin el eje axial del 
chorro. 

Se ha optado por la ultima de las dos posibihdades mencionadas, y el montaje 
experimental esta inspirado en los descriptos por Taub [1] y Wetsel [2]. El principio de 
funcionamiento consiste en detectar los cambios en el radio del chorro producidos por las 
perturbaciones, a traves de cambios en la sombra proyectada por el chorro al ser atravesado 
por un haz de laser. Este cambio en la sombra es aprovechado en un fotodiodo que genera una 
corriente en funcion de la intensidad de luz que recibe. Luego un tratamiento de las senales 
obtenidas permitira establecer el radio en funcion del tiempo para cada una de las posiciones 
elegidas y valores de tension aplicados en el electrodo. 

El liquido a utilizar en las experiencias sera agua destilada, que al estar exenta de 
particulas de tamaiio considerable, facilita que el chorro se forme en el orificio con el minimo 
de perturbaciones libres. 

El montaje experimental es dividido en las siguientes partes, a los fines de su 
descripcion: 

• Sistema de formacion del chorro. 
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• Dispositivo de desplazamiento espacial del chorro. 

• Sistema de perturbacion del chorro. 

• Sistema de electrificacion del chorro. 

• Sistema de medicion del radio del chorro y adquisicion de senales. 

Se detallan a continuacion cada una de las partes. Para el caso del sistema de medidas, 
se incluyen los calculos para la determinacion experimental de la escala entre los valores de 
tension medidos y las amplitudes de las perturbaciones en la superficie del chorro. 

Posteriormente, en el anexo de este capitulo se incluyen las especificaciones mas 
relevantes de algunos de los elementos que forman parte del montaje experimental. 

3.2 - SISTEMA DE FORMACION DEL CHORRO 

Para describir este sistema se comienza per detallar el sistema hidraulico que sirve para 
establecer el escurrimiento. Un esquema del mismo puede verse en la figura a continuacion. 






0 



0 



10 



RECIPIENTE NITROGENO 
VALVULA DE PASO 
REGULADOR DE RRESION 
VALVULA DE SEGURIDAD 



6 - BOCA DE CARGA Y PURGA DE AIRE 

7 - VALVULA DE AGUJA 



8 Fl OTAMF^RO 



9 - EILTRO D: YEMBRANA 



RECIRIENTE DE AGUA 



10 - BOOUILLA 



Figura 3.1 - Esquema simbolico del sistema hidraulico. 
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El liquido es almacenado en un recipiente a presion (5), construido en acero inoxidable 
AISI 304. El recipiente posee una capacidad de almacenamiento de aproximadamente un litro. 
Esta capacidad, para los caudales manejados, resulta en una autonomia de funcionamiento de 
aproximadamente una hora. 

Para la formacion del chorro se presuriza al liquido empleando N2, contenido en el 
recipiente apropiado (1), con sus correspondientes valvulas de paso (2) y de regulacion de 
presion (3), la cual incluye los manometros de alta y baja presion. La presion de operacion del 
sistema se fija en aproximadamente 50 kPa. Para evitar sobrepresiones que pongan en riesgo 
al sistema, entre el recipiente de liquido y el de N2 se coloca una valvula de seguridad (4), con 
una presion de apertura de 300 kPa, correspondiente con la presion admisible por el elemento 
mas debil (porta filtro de polipropileno). La carga del liquido al sistema, se realiza en forma 
manual desde la parte superior del recipiente, a traves de la boquilla (6). 

Para efectuar la medicion y regulacion del caudal volumetrico del chorro, 
posteriormente al recipiente se halla un medidor de caudal de area variable o flotametro (8), 
que tiene incorporada una valvula de aguja de accionamiento manual, la cual posee una muy 
alta sensibilidad (7). Este elemento va dispuesto sobre un soporte metalico que posee cuatro 
tomillos a los efectos de conseguir su nivelacion. 

Luego del caudalimetro, se encuentra un filtro de membrana de fibras de vidrio montado 
en un porta filtros de polipropileno (9). El filtro, unido al uso de agua destilada, elimina la 
posibilidad de taponamientos en el capilar en donde se ft)rma el chorro, y ruidos en las senales 
medidas que serian provocados por el pasaje de particulas a traves del orificio. Las 
membranas filtrantes empleadas poseen un diametro de 47 mm y un tamano nominal de poro 
de 1,5 |am. 
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El liquido es forzado a salir a traves de un orificio calibrado de 200 |u,m de diametro, 
formado por un capilar de acero inoxidable, adosado con pegamento al piezolelectrico que 
forma parte del sistema de excitacion del chorro, descripto mas adelante. El capilar es 
enrasado con las superficies del piezoelectrico usando una pulidora metalografica. El 
piezoelectrico va montado en una boquilla (10), construida tambien en acero inoxidable AISI 
304. 

En las figuras a continuacion se puede ver la construccion de la boquilla con detalle. 
Algunos de los aspectos que pueden notarse en las figuras seran considerados posteriormente 
en ocasion de describir el sistema de excitacion del chorro. 



4 




1 - CA^I.. AR A:FR0 INaX 5 - rUFRPG BQQJI I A 

2 - PIEZOELEFTRIFO £ - O'HINC 
3-0'RINCS 7 - A[CP_E MAKCUERA 
L - TUERCA 



Figura 3.2 - Corte de la boquilla y detalle de la union. 1 -Capilar de acero inoxidable. 
2-Piezoelectrico. 3-0'rings. 4-Tuerca de ajuste. 5-Cuerpo de boquilla. 6- O'ring. 

7-Acople a manguera. 
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Figura 3.3 - Perspectiva explotada de boquilla. 
Todos los componentes del sistema de formacion del chorro, se hallan unidos entre si 
mediante mangueras plasticas con malla de refuerzo de Nylon, que soportan presiones de 2 
MPa, y que son aseguradas a los acoples por abrazaderas metalicas. 

El sistema se completa con un dispositivo para recoger el agua, formado por un cano 
inclinado y acodado cuya descarga se situa en un recipiente, construidos en acero inoxidable. 

3.3 - DISPOSITIVO DE DESPLAZAMIENTO ESPACIAL DEL CHORRO 

Para poder seguir la evolucion de las perturbaciones en funcion del tiempo, en distintas 
posiciones segun el eje axial del chorro, se debe contar con un dispositivo que permita 
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cambiar el punto en donde el haz de laser intercepta al chorro. Los parrafos siguientes estan 
consagrados a la descripcion de tal sistema. 

La boquilla y su sistema de desplazamiento axial estan montados sobre una mesa de 
desplazamiento vertical sobre correderas, formando parte de un mecanismo provisto de una 
transmision por cadena, impulsado por un motor paso a paso. Estos elementos permiten 
realizar un ajuste fino de la altura del chorro, en correspondencia con el haz de laser, para 
eventualmente compensar la muy leve caida del chorro debida a la accion de la gravedad. Este 
ajuste solo puede ser necesario cuando se trabaja con muy bajos caudales y si se desea 
efectuar la medicion en puntos muy alejados del orificio. En la figura siguiente, pueden verse 
el motor paso a paso y su placa de montaje, la transmision compuesta de correa y cadena, la 
mesa de desplazamiento vertical y las barras que ofician de correderas. 




Figura 3.4 - Foto del dispositivo de desplazamiento vertical. 
El control del sentido de giro y del niimero de pasos a realizar del motor paso a paso, es 
efectuado a traves de un programa por medio del puerto paralelo de una PC. La salida de este 
control es tomada por los circuitos electronicos de potencia y control propios del motor. Estos 
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circuitos de potencia y control estan energizados por una fuente de tension independiente de 
la PC. 

El dispositivo de desplazamiento axial esta emplazado sobre la mesa de movimiento 
vertical, y esta constituido por una pequeiia plataforma en donde se sujeta la boquilla 
mediante una abrazadera. La plataforma puede desplazarse merced a un mecanismo de pinon 
y cremallera de accionamiento manual, al que se le ha adosado una regla graduada para medir 
la distancia entre la salida del chorro y la posicion en donde se mide su radio. 




Figura 3.5 - Sistema de desplazamiento axial. 




Figura 3.6 - Sistema de desplazamiento axial y escala de medicion de distancias. 
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Las figuras anteriores muestran dos vistas del dispositivo de desplazamiento axial. En la 
segunda de ellas puede apreciarse la regla graduada que permite la medicion de distancias. 

3.4 - SISTEMA DE EXCITACION DEL CHORRO 

El sistema de excitacion del chorro esta constituido por un ceramico piezoelectrico 
dispuesto sobre una placa circular de laton amarillo, en cuyo centre se aloja el capilar de acero 
inoxidable a traves del cual es establecido el escurrimiento. Para evitar la fuga de liquido 
alrededor del piezoelectrico, se emplean O 'rings que se situan en asientos conformados en la 
boquilla y en la tuerca que sostiene al conjunto (ver figuras 3.2 y 3.3). 

El piezoelectrico es excitado conectandolo a un generador de onda cuadrada, 
correspondiendo el borne positivo al ceramico y el negativo a la placa. El movimiento 
vibracional forzado del piezoelectrico, que tiene la direccion del eje axial del chorro, es el 
encargado de producir las perturbaciones en la superficie del mismo. 

3.5 - SISTEMA DE ELECTRIFICACION 

El chorro esta conectado a tierra por medio de una conexion derivada de la boquilla en 
donde es formado. La conexion de tierra puede ser apreciada en la figura 3.6. El campo 
electrico es aplicado al chorro a traves de un electrodo coaxial con el, constituido por un tubo 
de acrilico al cual se ha adosado una hoja de Mylar conductor (con deposito de oxido de 
estafio), en su interior. La hoja conductora se conecta a una fuente variable de alta tension, de 
corriente continua. Para que el haz de laser no sea afectado por su pasaje a traves del tubo y 
de la hoja conductora, se le han practicado ranuras a ambos elementos, obteniendose el pasaje 
libre del haz. Un esquema dimensional del electrodo se presenta en la figura a continuacion. 



99 



Montaje experimental 



HojQ o/e Mylar^ conduc Lor^ 





Figura 3.7 - Esquema del electrodo. 
La tension aplicada al electrodo es controlada mediante un multimetro desde la 
alimentacion con baja tension a la fuente, habiendose calibrado la misma previamente. 

La fuente de alta tension esta protegida por un gabinete plastico, y en su circuito de 
alimentacion se incluye un interruptor y un fusible para protegerla contra fallas. 



3.6 - SISTEMA DE MEDIDA Y ADQUISICION DE SENALES 

Todos los elementos del sistema optico estan montados sobre un banco construido a 
partir de dos bancadas de tomo, las cuales se apoyan sobre una mesa por medio de tomillos 
que permiten efectuar la regulacion de altura y la nivelacion. 

Sobre las bancadas se montan unas piezas de acero maquinadas en su parte inferior con 
la seccion transversal de la bancada; y en su parte superior, con una cola de milano que posee 
agujeros roscados en sentido vertical. Estas piezas se emplean para fijar o deslizar los soportes 
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de los postes que sujetan los elementos opticos. La fijacion de las piezas a la bancada, una vez 
posicionadas, se realiza a traves de tornillos. 




Figura 3.8 - Piezas de fijacion de los elementos opticos al banco. 
En primer termino se encuentra montado sobre el banco un laser que es vinculado al 
poste a traves de una mesa con desplazamiento angular segiin dos ejes y con carga por resorte, 
que posibilita el ajuste del haz de laser al eje optico del sistema. 

Luego del laser se dispone un sistema de dos lentes plano-convexas con el fm de 
aumentar el diametro del haz hasta un valor de aproximadamente 2 mm, tal que se cumpla una 
relacion aproximada de 1:10 entre el diametro del chorro sin perturbar y el diametro del haz. 
Cada una de las lentes esta fijada en un tambor, y ambos tambores se ajustan mediante acople 
roscado a un tubo de extension de 40 mm de longitud. Este tubo de extension permite regular 
la distancia entre las lentes y por ende el enfocado del sistema. El tubo de extension es 
vinculado al poste a traves de un anillo de fijacion. 

El haz de laser es luego desviado 90° para interceptar al chorro, empleando un espejo de 
22 X 22 mm, montado en una pequeiia base con desplazamiento angular segun dos ejes y con 
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carga por resorte. La base permite ajustar la direccion del haz reflejado, mediante tomillos de 
64 hilos por pulgada. En la otra bancada y enfrentado con el primero, otro espejo refleja 
nuevamente el haz 90° para alinearlo con el resto de los dispositivos opticos. Los espejos 
utilizados son del tipo aluminizados en la superficie mas cercana a la luz incidente, con una 
precision superficial correspondiente a % de onda a 632,8 nm, lo cual elimina las perdidas de 
energia y deformaciones en la imagen al evitar que el haz atraviese el vidrio. 

Una lente plano-convexa permite luego que la sombra proyectada por el chorro se 
magnifique en un factor de 10. Esto se debe a la necesidad de que la sombra tenga una altura 
de aproximadamente un tercio de la correspondiente a la ranura que se antepone al fotodiodo, 
y para que las lineas de difi-accion no pasen a traves de la mencionada ranura. La lente se 
monta en un anillo provisto de tres tomillos con extremos protegidos por Nylon, que facilitan 
su centrado. 

El sistema optico concluye en el fotodiodo, que es el encargado de generar las seiiales 
electricas de acuerdo a la intensidad de luz que recibe. 

El fotodiodo esta montado sobre una placa de acrilico, la cual esta sujetada a traves de 
cuatro tomillos, concentricos con un resorte cada uno de ellos, a una caja metalica formada 
por un cano cuadrado que posee una abertura vertical en uno de sus lados. El interior del cafio 
esta aislado de la luz ambiente por medio de tapas superior e inferior, y cintas adhesivas en su 
parte frontal. La presencia de los tornillos y resortes permite efectuar un correcto 
posicionamiento del fotodiodo respecto del haz de laser incidente, siendo este conducido al 
interior de la caja, atravesando una ranura de ancho ajustable por medio de tomillos adosados 
a la caja. La inclusion de esta ranura se justifica a partir de la necesidad de estudiar la 
evolucion en el tiempo del radio del chorro en una muy pequefia region del mismo, y a que 
esta Uegando a la caja un haz circular que abarca una region del chorro demasiado extensa. De 
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la caja se deriva la conexion electrica al sistema de adquisicion de datos y visualizacion de 
sefiales. Un detalle del conjunto descripto y un esquema del sistema optico pueden verse en 
las siguientes figuras. 




Figura 3.9 - Vista frontal y corte transversal del montaje del fotodiodo. 1-Caja metdlica. 
2-Placa de acrilico. 3-Fotodiodo. 4-Tornillos de ajuste ancho de ranura. 5-Placas de acero 
que forman la ranura. 6-Ajuste de posicion del fotodiodo (tornillos y resortes). 




Figura 3.10 - Representacion esquemdtica del sistema optico. 1 -Laser. 2 -Sistema de lentes 
conjugadas. 3-Espejos. 4-Lente de magnificacion. 5-Caja con fotodiodo y ranura. 
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La corriente generada por el fotodiodo en funcion de la intensidad de luz que incide 
sobre el, es luego amplificada a traves de un circuito electrometrico (que convierte la serial de 
entrada de corriente en una salida de tension), similar al utilizado por Desimone [3]. El 
diagrama simbolico de tal circuito es mostrado en la figura a continuacion. 
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Figura 3.11 - Esquema del circuito electrometrico. 
La seiial obtenida es luego amplificada 20 dB (o con una ganancia de 10), con el objeto 
de aumentar la relacion entre la serial de salida del circuito electrometrico y el ruido propio de 
la placa de adquisicion, de forma tal que la serial adquirida no se vea afectada, sobre todo 
cuando se amortigua. Esto evita inconvenientes a la hora de determinar el ancho de los pulsos. 

Un filtro pasa bandas, colocado posteriormente a este amplificador, permite eliminar los 
ruidos que alteran la serial. Las frecuencias de corte de pasa altos y pasa bajos se fijan 
respectivamente en 1 10 Hz y 4 kHz. 
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A la salida del filtro se dispone de un osciloscopio para el monitoreo continuo de las 
senales generadas, y una placa de adquisicion de datos colocada en una PC. 

3.6.1 - Estimacion de la escala tension - amplitud de perturbacion 

Para efectuar esta estimacion se requiere en primer termino relacionar la tension de 
salida del conjunto fotodiodo y circuito electrometrico, y luego tomar en consideracion que 
esa serial de tension se amplifica en un factor de 10 (20 dB) y que ademas la sombra 
proyectada por el chorro Uega al fotodiodo magnificada diez veces. 

La respuesta del sistema formado por el fotodiodo y el circuito electrometrico sigue la 
siguiente ley lineal [2]: 



Donde P es la potencia del haz de laser que incide sobre el fotodiodo, B es el factor de 
conversion del sistema entre intensidad de luz y tension, H es la altura de la ranura, Vd la 
tension residual, Rs el radio de la sombra del chorro sin perturbar y 5s la amplitud de la 
sombra de la perturbacion. 

La tension residual Vd se mide con la ranura cubierta, es decir con el fotodiodo a 
oscuras. El valor obtenido es Fc/ = - 0, 0033 V. 

Con la ranura plenamente iluminada con el haz de laser, puede determinarse el valor del 
producto B-PdL partir de la relacion siguiente: 




(3.1) 



V'=B-P + Vd 



(3.2) 
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Con el valor medido V'=3, 0471 V, y el anterior dado para Vd, resulta B P = 3,0504 ■ V . 

Midiendo la tension de salida cuando la ranura es iluminada por el haz ensombrecido 
per el chorro sin perturbar, puede determinarse la relacion entre la altura de la sombra del 
chorro y la altura de la ranura. El valor medido es V" = 1,8109 V, y la relacion buscada se 
halla a partir de: 



V"=B-P 



l-§yVd (3.3) 



Ds 

Resulta entonces un valor de — = 0,405 . Este valor concuerda en forma muy aceptable 

con el calculado a partir los valores fijados para la altura de la sombra del chorro y la de 
ranura, que son de 2 mm (recordar que el diametro del chorro sin perturbar es de 0,2 mm y la 
lente magnifica diez veces este valor), y 5 mm respectivamente. 

Luego, un valor de tension de 1 V medido en la placa de adquisicion de datos, se 
corresponde con una tension de 0,1 V en la salida del circuito electrometrico, y esta con un 
cambio en el radio de la sombra del chorro de 1,4153 mm, calculado con la expresion (3.1). 
Introduciendo la magnificacion de la lente, el cambio en el radio del chorro correspondiente 
sera de 0,14153 mm. 

Por lo tanto, los valores de tension en Volt que se obtienen a partir de la adquisicion de 
datos, deberan ser afectados por el factor 141,53 para ser Uevados a [om de variacion en el 
radio del chorro. 

Para finalizar con este apartado es menester efectuar la siguiente aclaracion. En realidad 
interesa conocer el radio de la sombra del chorro sin perturbar, la amplitud de la sombra de la 
perturbacion y la altura del rectangulo del haz de laser que inciden sobre el fotodiodo, luego 
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de atravesar la ranura, por la razon de que es el fotodiodo el elemento que genera las senales 
de corriente. No obstante, hay una relacion geometrica entre las magnitudes antedichas y las 
usadas en los calculos, correspondientes con la posicion y geometria de la ranura, relacion que 
esta dada por el angulo del cono que forman la sombra y el haz luego de pasar por la lente de 
magnificacion. Los calculos anteriores seran validos siempre y cuando el segmento del haz de 
laser que Uega al fotodiodo, luego de atravesar la ranura, este plenamente incluido dentro del 
area sensible del mismo. En este caso, tal condicion se cumple ya que luego de atravesar la 
ranura, la altura del segmento del haz de laser crece hasta 5,25 mm aproximadamente, siendo 
la altura del area sensible del fotodiodo, de unos 6,4 mm. 

3.6.2 - Precision del sistema de medida de perturbaciones 

En el dispositivo experimental descripto en [1], la precision del sistema aparece limitada 
por el ruido de las senales. Tal ruido esta ligado a las perturbaciones libres en la superficie del 
chorro y al ruido propio del sistema de visualizacion, y en ese trabajo se detalla que alcanza 
un valor del 0,5% del radio del chorro no perturbado. 

En el caso de este montaje experimental, dependiendo de la velocidad del flujo, se ha 
observado que el ruido oscila entre el 0,5 y el 1,2 % del radio del chorro sin perturbar, 
aproximadamente. Este ruido sera tenido en cuenta al momento de comprobar las tasas de 
crecimiento, segiin se vera. 
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3.7 - ESQUEMA Y VISTA DEL MONTAJE EN CONJUNTO 



A modo de conclusion de lo expresado a lo largo de este capitulo, se incluyen a 
continuacion una representacion esquematica general y una fotografia del montaje 
experimental. 




Figura 3.12 - Esquema general del montaje experimental. 1 -Laser. 2-Lentes conjugadas. 
3-Espejos. 4-Lente de magnificacion. 5-Ranura ajustable. 6-Fotodiodo. 7-Sistema recolector 
de agua. 8-Electrodo. 9-Boquilla, piezoelectrico y capilar. 10-Mesas de desplazamiento del 
chorro. 1 1-Motor paso a paso y transmisidn. 12-Control motor. 13-Fuente alimentacidn 
motor. 14-PC control motor. 15-Fuente alta tension. 16-Excitaci6n a fuente alta tension. 
17-Filtro de agua. 18-Flotdmetro y vdlvula de aguja. 19-Recipiente de agua. 20-Vdlvula de 
seguridad. 21-Circuito electrometrico. 22-Amplificador. 23-Filtro pasa bandas. 
24-Osciloscopio. 25-PC con placa de adquisicion de senales. 26-Generador de ondas. 
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ANEXQ 3.1 - ESPECIFICACIONES TECNICAS DE EQUIPQS 
Caudalimetro 

• Marca: GILMONT INSTRUMENTS. 

• Modelo: ACCUCALGF-4542-1230 

• Precision: ±1% 

• Material en contacto con el fluido: acero inoxidable. 

• Presion maxima de entrada: 200 psig. 

• Temperatura maxima del fluido: 250 °F. 

• Temperatura maxima ambiente: 200 °F. 
Generador de ondas 

• Marca: GOOT 

• Modelo: 1001 

• Rango de frecuencias: 4 - 200 kHz. 

• Rango de tension: 0 - 5 V (eficaz), para onda sinusoidal. 

0 - 15 V (pico a pico), para onda cuadrada. 
Fuente de alta tension 

• Marca: BRANDENBURG 

• Modelo: Serie 569 

• Rango de tension de entrada: 5-15 Vcc. 
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• Rango de tension de salida: 2-6 kV. 

• Corriente de salida a tension maxima: 0.5 mA. 

• Ruido de salida: 0.2% 

• Potencia: 3 W 

• Rendimiento: 70%. 
Laser 

• Marca: UNIPHASE 

• Tipo: Helio - Neon 

• Clase: Ilia 

• Potencia minima salida: 0,8 mW TEMoo 

• Longitud de onda: 632,8 nm 

• Diametro: 0,48 mm 

• Divergencia: 1,7 mrad 

• Ruido maximo (valor eficaz): 1% (30 Hz - 10 MHz) 
Fotodiodo 

• Para infrarrojos, de silicio, de alta velocidad y area media, en un encapsulado 
hermeticamente sellado de tipo T08. 

• Marca: IPL 

• Respuesta espectral de pico: 800 nm. 

• Longitud de onda: 350 a 1 100 nm 
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• Responsividad de pico: 0,5 A/W 

• Area activa: 41,3 mm (area sensible cuadrada) 

• Corriente residual: 4 nA (Vr = 1 V) 

• Capacitancia: 325 pF (Vr = 0 V) 

• Tension inversa de ruptura: 60 V 

• Tiempo de respuesta: 25 ns (Vr = 10 V, Rl = 100 Q) 

• Temperatura de funcionamiento: -40 a +70 "C 
Amplificador 

• Marca: HEWLETT - PACKARD 

• Modelo: 465 A 

• Ganancia de tension: 20 dB (xlO) o 40 dB (xlOO), circuito abierto. 

• Precision ganancia: ± 0,1 dB (± 1%) a 1000 Hz 

• Respuesta en frecuencia: ±0,1 dB, 100 Hz a 50 kHz; <2 dB, 5-100 Hz y 50 kHz-1 
MHZ 

• Salida: > 10 V eficaz, circuito abierto; > 5 V eficaz en 50 Q 

• Distorsion: 1%, 10 Hz a 100 kHz; < 2%, 5 Hz a 10 Hz y 100 kHz a 1 MHz 

• Impedancia de entrada: 10 MQ, < 20 pF 

• Impedancia de salida: 50 Q 

• Ruido: < 25 ^iV eficaz referido a entrada con 1 MQ resistencia de fiiente. 
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Filtro pasa bandas 

• Marca: KROHN - HITE 

• Modelo: 3100R 

• Rango: 10 Hz - 1 MHz, regulacion independiente de frecuencia de corte alta y baja. 

• Ganancia: 0 dB 

• Tasa de atenuacion fuera de la banda: 24 dB/octava 

• Tasa de atenuacion maxima: 80 dB 

• Ruido maximo: 100 [iV 

• Tension maxima entrada: 3 V (eficaz) 
Placa de adquisicion de datos 

• Marca: ADVANTECH CO., LTD. 

• Modelo: PCL 1800 

• Velocidad maxima de muestreo: 200 kHz. 

• Entrada analogica: Canales: 16 

Resolucion: 12 bit 

Tiempo de conversion: 3 \is 

Tension de entrada: ± 10 V (rango maximo) 

Disparo: por software (incluido con la placa) 

Precision: 0.01% FSR + 1 LSB. 
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4 - ESTUDIO EXPERIMENTAL 
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4.1 - EVTRODUCCION 

Este capitulo tiene como objeto la presentacion de los resultados experimentales 
obtenidos mediante el montaje anteriormente descripto, y la comparacion de los mismos 
frente a las predicciones teoricas. 

Se comienza por establecer los limites sobre algunos de los parametros experimentales, 
los cuales permiten justificar la eleccion de sus valores. Luego se presenta una serial 
caracteristica de las obtenidas con el sistema de adquisicion, para posteriormente, y con apoyo 
en ella, describir la manera en la cual se les ha dado tratamiento para la determinacion de las 
tasas de crecimiento y de las extensiones del paquete de ondas. Finahnente, se presenta un 
conjunto de seiiales caracteristicas y se muestran los resultados experimentales obtenidos para 
las tasas de crecimiento y los limites del pulso, realizando a su vez una comparacion de estos 
y de las caracteristicas de las seiiales con las predicciones teoricas, lo cual permite extraer 
conclusiones acerca de la validez del modelo presentado previamente. 

4.2 - LIMITES EXPERIMENTALES 

En este apartado se discuten las limitaciones del estudio experimental asociadas con la 
velocidad del chorro, las posiciones segun la coordenada axial en donde se efectuan las 
medidas, la excitacion acustica y la excitacion electrica. Estas limitaciones estan relacionadas 
con el montaje experimental algunas, y con la dinamica del escurrimiento otras. Las 
consideraciones aqui presentadas son las que justifican los valores de dichos parametros 
experimentales, adoptados en este estudio. 
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4.2.1 - Limites de velocidad del chorro 

En primera instancia, este estudio esta limitado al campo de inestabilidades convectivas, 
ya que cuando la inestabilidad es absoluta, el problema de la senal no es posible debido a que 
no puede establecerse la excitacion con un patron de ondas. Teoricamente, el efecto del 
campo electrico sobre el escurrimiento es mas notorio en cercanias de la transicion entre uno 
y otro tipo de inestabilidad, esto es para bajos valores de We, lo cual se corresponde para el 
liquido y diametro del orificio elegidos, con velocidades bajas del flujo. 

La primera limitacion a la velocidad del chorro aparece asociada con efectos 
gravitatorios. Cuando la relacion entre las fuerzas de inercia y las gravitatorias es baja, el 
chorro comienza a caer notoriamente en cercanias del orificio perdiendo la coaxialidad con el 
electrodo, lo cual promueve dificultades a la hora de efectuar las medidas y es una alteracion a 
la geometria estudiada en forma teorica. Ademas, el hecho de que el chorro sea afectado por 
efectos gravitatorios, constituye un alejamiento de las hipotesis que sustentan el estudio 
teorico. Debe recordarse que el estudio teorico presentado en el capitulo 2 es uno local (el 
perfil de velocidad constante es invariante con la coordenada axial), y que el analisis teorico 
de un chorro sometido a fuerzas gravitatorias requiere de los metodos globales, ya que todas 
las estaciones segiin la coordenada z no tendrian el mismo perfil de velocidad. El valor mas 
bajo de velocidad empleado es de Uo=7,2 m/s, que se corresponde con We=72. 

Cuando la velocidad del chorro es alta, el mismo se vuelve en teoria menos inestable y 
los efectos de la excitacion electrica sobre la propagacion del pulso son en teoria menos 
notorios. El valor mas alto de velocidad empleado es de Uo= 10,6 m/s, que resulta en We=156. 
Se ha adoptado este valor, de forma tal que no sea necesario aumentar en demasia la tension 
de excitacion del piezoelectrico para mantener una relacion aceptable entre senal y ruido (lo 
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cual puede Uevar a la destruccion del actuador o bien producir senales anormales), y para no 
afectar el funcionamiento del piezoelectrico a causa de someterlo a una presion elevada por 
parte del liquido. 

El estudio se completa con un tercer valor de velocidad, Uo=9,9 m/s, para el cual se 

tiene We=136. 

4.2.2 - Limitaciones a las posiciones de medicion segun la coordenada axial 

Existe un limite inferior para la posicion axial a partir de la cual pueden realizarse las 
medidas, que esta ligado a la construccion de la boquilla y a la dimension del haz de laser en 
su interseccion con el chorro. Este limite es de z/=5 mm, en donde no hay interaccion entre el 
laser y la pieza roscada que sostiene al piezoelectrico contra la boquilla. 

La teoria desarroUada para el fenomeno en estudio se basa en una aproximacion lineal 
de la superficie del chorro perturbado. Se sabe que cercano a la zona de ruptura, en donde la 
amplitud de las perturbaciones es importante en relacion al radio, los efectos no lineales 
toman relevancia. Aqui se toma un limite superior para la posicion de medicion de Ls=ll mm, 
para la cual la amplitud maxima de la perturbacion alcanza, en el caso mas desfavorable, un 
valor comprendido entre el 10 y el 15 % del radio del chorro sin perturbar. 

Las medidas son efectuadas entre los limites antedichos con un salto zlz de 1 o 2 mm 
entre cada una de las posiciones. Esto proporciona varios puntos diferentes de observacion, 
que se consideran suficientes para analizar la evolucion del paquete de ondas. 
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4.2.3 - Limites de la frecuencia de excitacion acustica 

El montaje experimental pennite obtener la respuesta del chorro en funcion del tiempo 
para cada posicion fija segun la coordenada axial. Un grafico tridimensional del radio en 
funcion de la posicion axial y del tiempo puede ser construido a partir de los datos obtenidos 
en las posiciones elegidas. La figura 4. 1 es una representacion esquematica de la propagacion 
de las perturbaciones en el dominio z-t, en la cual se muestra la senal obtenida para una 
posicion fija, ademas de la reconstruccion de las perturbaciones para dos tiempos diferentes. 




^ Coordenada 
L axial z 



Figura 4.1 — Esquema de la propagacion de una perturbacion pulsante. 
Las lineas identificadas por z/tj y z/tn indican los rayos de comienzo y final del paquete 
de ondas. Por fiiera de la region delimitada por estos rayos las perturbaciones son 
amortiguadas, y son amplificadas dentro de ella. Las mediciones deben ser efectuadas en 
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posiciones donde no exista interaccion entre dos pulsos sucesivos, de forma tal que el ancho 
del paquete de ondas pueda ser determinado. 

Observando la figura puede deducirse facilmente que para satisfacer la mencionada 
condicion, para un chorro de una longitud L (o para la mayor distancia entre la fuente de 
perturbacion y la posicion de observacion), el tiempo transcurrido entre la excitacion de dos 
pulsos sucesivos debera ser mayor que el intervalo At^^^ dado por [1]: 



Por ejemplo, empleando las predicciones teoricas ya descriptas, para un chorro de radio 
a=100 /Mn, We=72, y una distancia L=ll mm, vale At^^^ =0,62-ms. La frecuencia de la 

serial cuadrada impuesta al piezoelectrico debera ser entonces menor a f=1575 Hz, para este 
caso. 

4.2.4 - Limite de la excitacion electrica. Descarga corona. 

La descarga corona se caracteriza por la produccion de iones en un fluido dielectrico (en 
este caso aire), y se produce cuando el valor del campo electrico es lo suficientemente intense. 
Esta nube de iones migra en el fluido dielectrico, bajo el efecto de fuerzas coulombianas, 
alcanzando una distribucion estable en el espacio luego de un cierto tiempo. Este fenomeno 
hace que el campo electrico pierda su configuracion inicial, y por lo tanto la superficie del 
chorro no resulta excitada por el mismo en la manera que se preve. La aparicion de la 
descarga corona establece entonces el limite para la excitacion electrica. 




(4.1) 
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Para una geometria de dos electrodos cilindricos, coaxiales e infinitos, la expresion 
experimental que permite calcular el valor del campo electrico a partir del cual comienza la 
descarga corona, para aire a presion y temperatura ambiente, se escribe [2]: 



E,„Xvlm)=3,l-10 



6 



^ J _^ 0,308 ^ 



(4.2) 



donde a es el radio del chorro. 

En este caso, con el radio del chorro fijado en a=0,01cm, se obtiene 

^cor = ^^'^ ' Vjin . A partir de la expresion que relaciona el valor del campo electrico con 

la diferencia de potencial, para una geometria cilindrica infmita, se encuentra el limite de la 
tension de excitacion: 



V^„^ = - a- In 



max cor 



= 12,6- 10" -10 ' -In 



^350^ 



V^7380-V 



en la cual b es el radio interior del electrodo. 

Debido a que no se trabaja con una geometria infinita, y a que posiblemente existan 
pequenos alejamientos de la geometria cilindrica coaxial, se toma un margen de seguridad y 
se fija el valor maximo de tension en Vmax=6 kV. El valor minimo de tension empleado es 
Vmin=2 kV, y el estudio se completa con tensiones fijadas en 3, 4 y 5 kV. 

4.3 - TRATAMIENTO DE SENALES 

En esta seccion se detallan los procedimientos usados para la determinacion de la 
extension del paquete de ondas y de las amplitudes de las perturbaciones, a partir de las 
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senales medidas, y el tratamiento seguido para establecer la tasas de crecimiento empleando 
dichas amplitudes. 

Para cada posicion axial, valor de campo electrico y de velocidad, se han obtenido tres 
muestras con el sistema de adquisicion. Cada una de esas muestras esta compuesta por varios 
pulsos, de los cuales se han seleccionado los que presentaron un comportamiento mas claro y 
han sido menos afectados por el ruido, respecto del resto. 

En adelanto a la presentacion de los resultados experimentales, se incluye aqui en la 
figura 4.2, una de las senales tipicas obtenidas mediante el sistema de adquisicion, con el 
objeto de facilitar la descripcion de estos metodos. 




-6- 



0.015 0.020 0.025 0.030 

Tiempo (ms) 

Figura 4.2 - Serial caracteristica. Frecuencia de excitacion acustica: 25 Hz aprox. 
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La senal mostrada ha sido afectada por la escala entre tension y amplitud calculada en el 
capitulo anterior. Cualquier comentario que pudiera tiacerse sobre la misma, en relacion con 
lo esperado en funcion de las previsiones teoricas, sera dejado de lado por el momento, hasta 
la discusion de resultados. 

En la figura puede observarse la serial tal cual es muestreada con el sistema de 
adquisicion (recordar que se han fijado las frecuencias del filtro pasa bandas en 110 Hz y 4 
kHz), y puede observarse el filtrado digital de la serial que elimina el ruido propio de la placa 
de adquisicion (este filtrado reconstruye la serial, a partir de los datos tornados cada 15 puntos 
de la sefial original). En adelante, todas las figuras de sefiales presentadas corresponden a la 
salida del filtro digital. 

4.3.1 - Determinacion de las tasas de crecimiento 

De acuerdo a las predicciones teoricas, las amplitudes de las perturbaciones 
correspondientes a los distintos rayos z/t, crecen de forma exponencial con una tasa de 
crecimiento temporal dada por la ecuacion (2.9). La amplitud de las perturbaciones pueden 
entonces ser descriptas mediante la siguiente ecuacion, con r/o la amplitud inicial y cr la tasa 
de crecimiento correspondiente con un rayo z/t dado: 



Tj = JJo-e 



at 



(4.3) 



Considerando la tasa de crecimiento adimensional y el tiempo expresados como: 



z 



a 



t = 




a = cr- 



U, 



0 



la ecuacion anterior puede escribirse de la siguiente forma: 
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* z 

ri = ri,-e " (4.4) 

Las tasas de crecimiento, en forma adimensionalizada, se pueden obtener entonces a 
partir de un analisis de las senales caracteristicas para diferentes posiciones axiales . 

El estudio de las tasas de crecimiento ha sido efectuado para tres rayos z/t diferentes, 
correspondientes a tres maximos distintos de la senal, a saber: el de mayor amplitud positiva, 
el del ciclo que le sigue, y el de mayor amplitud negativa (A, B y C respectivamente, en la 
figura 4.2). Se determino de esta manera que el pico de mayor tasa de crecimiento es el que 
alcanza la maxima amplitud positiva (indicado en la figura 4.2 con A). Este pico es el que se 
corresponde con la mayor tasa de crecimiento y en consecuencia, teoricamente, con el rayo 
z/t=l. 

La determinacion de las tasas de crecimiento mediante este procedimiento, presenta dos 
aspectos que provocan una incertidumbre en los resultados obtenidos. Un primer aspecto esta 
asociado con el ruido de la serial. La obtencion de los valores de las amplitudes aludidas 
anteriormente se puede realizar directamente mediante lectura desde los graficos de las 
seriales. Las valores asi determinados estaran acompanados por una banda de error, que se fija 
en relacion al ruido de la senal, tambien leido desde los graficos. El segundo aspecto esta 
vinculado a la indeterminacion que existe en la medida de la posicion axial. Las bandas de 
error para la distancia axial al orificio se establecen por la apreciacion que se tiene con la 
regla graduada que se emplea para medirla. 

4.3.2 - Determinacion de la extension del paquete de ondas 

La presencia de ruido en las seriales, hace que no puedan establecerse con precision los 
limites del paquete de ondas en forma directa midiendo sobre los graficos de las senales. 
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sobre todo en la "cola" del pulso, donde la amplitud se atenua lentamente confundiendose con 
el ruido. Conocida entonces esta limitacion, se ha optado per determinar la extension del 
paquete de ondas entre el comienzo del mismo y el pico de mayor amplitud, midiendo sobre 
los graficos de las seiiales, ya que el comienzo del pulso resulta visible con relativa facilidad. 
Este procedimiento se repite para las tres senales que corresponden a una posicion, velocidad 
y excitacion electrica dadas, y el valor a tener en cuenta en cada caso resulta de promediar los 
tres obtenidos. 

Esos resultados pueden ser comparados con las previsiones teoricas, considerando a los 
rayos z/ti y z/t=l para el calculo de la extension del paquete de ondas segun la ecuacion 
(2.14). 

Un procedimiento altemativo al anterior y tal vez mas interesante, ya que permitiria 
establecer la extension del paquete de ondas completo, ha sido probado sin exito para la 
totalidad de las senales (se ha logrado que funcione solo para el valor mas alto de We). El 
mismo consiste en ajustar la serial mediante una funcion y a partir de ella, con un criterio 
dado, estimar el valor buscado del ancho del pulso. 

Observando la forma de la senal puede verse que se asemeja a la respuesta impulso de 
un sistema de segundo orden, segun la defmicion del mismo en la teoria de control clasica [3]. 
La respuesta de tal tipo de sistema en el dominio del tiempo es bien conocida, y aqui se ha 
empleado esa funcion para ajustar las sefiales, modificada y con el agregado de algunas 
constantes que tienden a una mejora en los resultados. 

La funcion de ajuste tiene la siguiente expresion: 

y(t) = yg+A-(l-e /-e ■sen((o-t + ^) (4.5) 
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en dondej^o. A to, xi, P,X2, co y ^son los parametros de ajuste. 

Para establecer el ancho del pulso puede trabajarse solamente con la envolvente de la 
funcion, que esta dada por el producto de las exponenciales y la constante A. EI criterio 
adoptado es que el pulso se ha amortiguado cuando la envolvente cae por debajo del 2% de su 
valor maximo. 

4.4 - RESULT ADOS Y DISCUSION 

Este apartado esta consagrado a la presentacion de las senales caracteristicas y los 
resultados obtenidos mediante los procedimientos detallados con anterioridad. En cada caso, 
la discusion y comparacion con la previsiones teoricas del capitulo 2 son efectuadas, y para 
ello se ha prestado atencion a las consideraciones del trabajo de Artana et al [4]. 

Los resultados caracteristicos para cada valor de campo electrico son presentados aqui 
en funcion del numero de Rayleigh electrico (el cual establece una relacion entre las fuerzas 
electricas y las asociadas a la tension superficial), definido mediante: 

Ra^=Eu^-We = p-E,J -"^ (4.6) 

7 

4.4.1 - Senales caracteristicas 

En la figura a continuacion se presentan resultados tipicos de la forma que adquiere la 
superficie del chorro perturbado, representada por las seriales adquiridas y luego filtradas 
digitalmente. Estas seiiales representativas han sido obtenidas para una velocidad del chorro 
de 7,2 m/s, a distancias del orificio de 7 y 10 mm, con una frecuencia de excitacion al 
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piezoelectrico de aproximadamente 25 Hz, y en el caso del chorro electrificado, con una 
tension aplicada al electrodo &q 5 kV. 



E 

■= 0 



We=72, Rae=0, 7=1 mm 




0.020 Q025 

lienipo (ms) 
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0.055 0.060 0.065 
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Figura 4.3 - Amplitud de la perturbacion sobre la superficie del chorro en funcion del 

tiempo. 

4.4.1,1 - Efecto del movimiento del piezoelectrico 



Puede apreciarse en la figura que el pulso siguiente al mostrado en forma completa, 
toma una amplitud de signo diferente en el inicio. Tal comportamiento esta ligado al hecho de 
que el movimiento vibracional del piezoelectrico, que posee la direccion del eje axial del 
chorro, tiene el mismo sentido que la velocidad del escurrimiento en uno de sus semi ciclos de 
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oscilacion, y contrario en el otro semi ciclo. Las senales caracteristicas en todos los casos 
analizados para establecer las tasas de crecimiento y el ancho del paquete de ondas, estan 
representadas por un pulso cuya amplitud inicial crece negativamente. No obstante ello, 
ninguna diferencia, en cuanto a los resultados, apareceria asociada al empleo de los pulsos 
cuya amplitud inicial crece positivamente, ya que se observa un comportamiento simetrico de 
ambos pulsos, respecto de la situacion no perturbada (la cual esta representada por una senal 
continua, formada por una linea recta paralela al eje de abscisas y ligeramente desplazada del 
origen del eje de ordenadas). 

4.4.1.2 - Forma de las senales caracteristicas 

Un hecho de gran importancia merece ser destacado: existe una asimetria respecto de la 
amplitud maxima de las perturbaciones, la cual como se ha puntualizado corresponde a z/t=l, 
teniendo en cuenta el tiempo que la separa del comienzo y final del paquete de ondas. Las 
predicciones teoricas concuerdan con esta observacion. En efecto, la simetria de los rayos z/ti 
y z/tn, respecto del rayo z/t=l, arroja como resultado que el tiempo que separa el comienzo 
del paquete de ondas, del pico que posee la mayor tasa de crecimiento, es menor que el 
correspondiente a la separacion de tal pico con la "cola" del paquete de ondas. Esta prediccion 
puede justificarse a partir de la expresion (2.14) cuando se reemplazan z/ti y z/tn por los 
valores apropiados (z/f/^ y z/tu'^ son una medida de la tangente con la cual las amplitudes 
crecen y decrecen, respectivamente, en la cabeza y cola del paquete de ondas). 

4.4.1.3 - Amplitudes de las senales caracteristicas 

Del analisis de la figura 4.3, puede destacarse que en los casos del chorro electrificado 
se observa un crecimiento en la amplitud maxima de la perturbacion, y que no se puede 
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apreciar (a simple vista) que el paquete de ondas resulte expandido con tal excitacion, dada la 
presencia del ruido. Puede observarse ademas que al aumentar la distancia al orificio, se 
produce un crecimiento en la amplitud de las perturbaciones junto con una expansion del 
paquete de ondas (si bien para esto ultimo no puede a simple vista determinarse claramente en 
que magnitud lo hace, se tiene una percepcion del fenomeno). Las previsiones teoricas 
concuerdan, al menos por el momento cualitativamente, con estos resultados. 

4.4.1,4 - Espectro de frecuencias de las senales caracteristicas 

En la siguiente figura se muestra el espectro de frecuencias de la serial mostrada en la 
flgura 2.2, obtenida mediante una transformada rapida de Fourier. 

Frequencia (Hz) 

0 500 1000 1500 2000 2500 3000 3500 4000 
0.5-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.4- 

0.3- 



0.2- 




0 500 1000 1500 2000 2500 3000 3500 4000 

Frequencia (Hz) 



Figura 4.4 - Espectro de frecuencias caracteristico. 
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Puede observarse una frecuencia dominante con un valor de aproximadamente 1350 Hz. 
Esta frecuencia, es una caracteristica de todas las senales muestreadas, con una variacion de 
+ 10 Hz, para todos los valores de velocidad y campo electrico empleados. Los cambios 
previstos para la extension del paquete de ondas, deben manifestarse entonces mediante la 
aparicion de una mayor cantidad de ciclos o semi ciclos en las senales, con anterioridad y 
posterioridad al pico de maxima amplitud. 

4.4.2 - Tasas de crecimiento maximas 

En las siguientes figuras se presentan, a modo de ejemplo, los graficos de la amplitud 
maxima de las perturbaciones en funcion de la distancia al orificio, en forma adimensional. 
Cada uno de estos graficos corresponde a un valor dado de velocidad y campo electrico. En 
los mismos, pueden ser observadas las barras de error asociadas al ruido en el caso de la 
amplitud, y a la precision de la medida de distancias en el caso de la coordenada axial. Las 
curvas en linea punteada corresponden al ajuste de los datos mediante una funcion como la 
mostrada en (4.4), donde se han dejado como parametros libres en el ajuste por cuadrados 
minimos tanto a la tasa de crecimiento como a la amplitud inicial. En cambio en linea 
continua se representa el ajuste con la misma funcion, pero fijando el valor de la tasa de 
crecimiento segun lo previsto en forma teorica, de forma tal que el unico parametro de ajuste 
es la amplitud inicial de la perturbacion. 
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Figura 4.5 — Amplitud maxima de la perturbacion en funcion de z para We-72 y Rae-0. 
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Figura 4.6 - Amplitud maxima de la perturbacion en funcion de z para We=136 y Rae=0. 
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We=156, Rae=0,58 

■ Amplitud medida 
Ajuste tasa teorica 




40 50 60 70 80 90 100 110 120 

z adimensional 

Figura 4.7 - Amplitud maxima de la perturbacion en funcion de z para We=156y Rae=0,58. 

La comparacion directa entre las amplitudes para diferentes valores de velocidad, no es 
posible ya que la tension de excitacion del piezoelectrico ha sido cambiada a efectos de 
reducir la importancia del ruido. El resto de los graficos empleados para la determinacion de 
las tasas de crecimiento son mostrados en el anexo 4.1. 

El crecimiento de la amplitud con respecto al aumento de la distancia al orificio, puede 
observarse en la gran mayoria de los casos. Cuando esto no acontece, probablemente se 
relacione con un cambio en la velocidad del escurrimiento, ya que la regulacion de caudal ha 
sido la situacion mas problematica al momento de realizar las mediciones. La informacion 
obtenida de todos estos graficos, en cuanto a la tasa de crecimiento en ftincion de la excitacion 
electrica, se resume en las figuras que se presentan a continuacion. En las mismas, las barras 
de error para las tasas de crecimiento representan el error obtenido en el ajuste. 
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Figura 4.8 - Tasa de crecimiento maxima en funcion de Rae, para We=72. 
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Figura 4.9 - Tasa de crecimiento maxima en funcion de Rae, para We=136. 
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Figura 4.10 - Tasa de crecimiento maxima en funcion de Rae, para We=156. 
En tanto que para We=156 y We=136 se observa una dependencia monotonica de la 
tasa de crecimiento con respecto a la excitacion electrica, no puede decirse lo mismo para 
We=72. Posiblemente el comportamiento observado para la velocidad mas baja empleada este 
asociado con el efecto de la gravedad, ya que el chorro en este caso evidencia una caida en 
posiciones mas cercanas al orificio, perdiendo la coaxialidad con el electrodo, y resultando de 
esta manera excitado de forma diferente a la prevista. No debe olvidarse ademas que, dadas 
las hipotesis del estudio teorico, el efecto de la gravedad es ignorado. 

Dejando de lado entonces los resultados experimentales obtenidos para la velocidad mas 
baja, puede decirse que la teoria concuerda cualitativamente con las observaciones 
experimentales. Valores mas bajos de las tasas de crecimiento maximas han sido observados, 
y ello tal vez este asociado a que el modelo teorico supone un chorro no viscoso, de tal forma 
que el efecto amortiguante de las fuerzas viscosas no es tenido en cuenta. 
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Otro hecho a destacar es que si bien para los valores mas altos del numero de Rayleigh 
se manifiesta un aumento en la tasa de crecimiento, el mismo resulta menor que el previsto. 
Una posible explicacion de este fenomeno es que, al estar la superficie del chorro perturbada 
por los pulsos, la configuracion del campo electrico es alterada, pudiendo iniciarse la descarga 
corona con un valor de tension menor que el previsto mediante la ecuacion (4.2). Otra posible 
causa de este fenomeno, aunque menos probable, es que la evolucion de las perturbaciones 
puede ser afectada por efectos no lineales, cuando la excitacion electrica alcanza un valor 
elevado. Recordando que el modelo teorico presentado en el capitulo 2 se basa en una 
descripcion matematica lineal de la superficie del chorro perturbado, se justificaria asi tal 
discrepancia. 

En relacion al comportamiento de la tasa de crecimiento maxima respecto de la 
velocidad del escurrimiento, los resultados experimentales muestran que las predicciones 
teoricas dan una tendencia acertada, en cuanto a que el escurrimiento es mas estable a medida 
que su velocidad aumenta. 

En resumen de lo expuesto acerca del efecto de la excitacion electrica sobre las tasas de 
crecimiento maximas, puede decirse que la teoria concuerda cualitativamente con los 
resultados experimentales que pueden considerarse como validos. Valores mas pequeiios de 
las tasas de crecimiento, y un efecto menos pronunciado que el previsto sobre las mismas 
cuando la excitacion electrica es fuerte, han sido obtenidos experimentalmente. Mediante este 
analisis se pone de manifiesto la necesidad de refinar la teoria, tomando en cuenta la 
viscosidad del liquido y tal vez aproximando la superficie del chorro perturbado mediante un 
modelo no lineal. En cuanto al efecto de la velocidad sobre la estabilidad del escurrimiento, la 
teoria concuerda cualitativamente con las observaciones experimentales. El aumento en la 



135 



Estudio experimental 



velocidad hace que el escurrimiento se comporte de manera mas estable. Las discrepancias 
encontradas en este caso tambien pueden estar asociadas a efectos viscosos. 

4.4.3 - Limites del paquete de ondas 

4.4.3.1 - Medicion direct a sobre elsrdfico 

Los resultados de las mediciones de la extension del paquete de ondas, entre el 
comienzo del mismo y la amplitud maxima, en funcion de la excitacion electrica, son 
mostrados en la tabla a continuacion y comparados con los previstos teoricamente (estos 
ultimos calculados con la expresion (2.14) empleando un valor de z/tu=l). Los datos que se 
muestran en esta tabla corresponden a We=72, y a z=5 mm. 



Rae 


Semiancho (ms) 


Teorico 


Medido 


0.00 


0,1731 


1,08 


0.14 


0,1744 


1,09 


0.32 


0,1769 


1,07 


0.58 


0,1823 


1,08 


0.90 


0,1917 


1,07 



Tabla 4.1 - Semiancho de los pulsos para We=72 y z=5 mm. 
Se observa que el semi ancho medido sobre las senales es un orden de magnitud mayor 
que el previsto por la teoria, y puede decirse ademas que alcanza un valor de 1,08 + 0,01 ms. 
El mismo valor del semi ancho ha sido observado para todas las posiciones, velocidades, y 
campos electricos, poniendo de manifiesto las dificultades que acarrea el ruido en la 
determinacion del mismo (si un cambio en el semi ancho del pulso realmente acontece en 
funcion de la excitacion electrica al chorro, es altamente probable que este comprendido 
dentro del ruido de la seiial). No se incluye aqui la comparacion entre las previsiones teoricas 
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y los resultados obtenidos para los otros valores de We testeados, ya que en esos casos 
solamente a partir de Rae»l,5 se preven cambios mmimos en la extension del paquete de 
ondas, y el estudio experimental esta limitado a Rae=l,3. 

4,4,3.2 - Medicion a partir de la funcion de afuste 
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Figura 4.11 - Resultado caracteristico del ajuste de senales. 
En la figura 4.1 1 se muestra a titulo de ejemplo una senal ajustada mediante la ecuacion 
(4.5), la cual no ha sido afectada por la escala entre tension y amplitud en |u,m. La misma 
corresponde a We=156, z=5 mm, sin excitacion electrica. Si bien el uso de este ajuste parece 
mas adecuado para Uevar a cabo un analisis del limite del paquete de ondas, lamentablemente, 
del conjunto de sefiales obtenidas, solamente para el We indicado el mismo pudo ser Uevado a 
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cabo satisfactoriamente. Para el resto de los We en cuestion, el ajuste no pudo realizarse ya 
que no se ha logrado la convergencia de los distintos parametros de la funcion de ajuste. 

A partir de los ajustes de las senales muestreadas para las posiciones axiales analizadas, 
el siguiente grafico ha sido construido, en donde se muestra la evolucion del ancho completo 
de paquete de ondas en funcion de z y se lo compara con las predicciones teoricas. 
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Figura 4.12 - Extension del paquete de ondas en funcion de la distancia al orificio. 

Puede observarse en esta figura que si bien el ancho de los pulsos es un orden de 
magnitud mayor que el correspondiente al previsto en forma teorica, existe una concordancia 
al menos cualitativa en cuanto a su aumento, en la medida que los pulsos se alejan del orificio. 
La extension de los pulsos, obtenida experimentalmente, muestra un crecimiento menor al 
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esperado, lo cual se acentua para posiciones mas alejadas del orificio. Este comportamiento 
tal vez este ligado a que en el ajuste, por la forma en que ha sido efectuado, no es tenida en 
cuenta la expansion de la senal en su "cabeza" (nuevamente aqui se manifiestan las 
dificultades ocasionadas por el ruido, que no permite establecer claramente el comienzo del 
paquete de ondas). No obstante, la influencia de los efectos viscosos y de la aparicion de 
efectos no lineales para posiciones alejadas del orificio, no deben ser descartadas. 

En teoria, para este valor de We no existen cambios para el ancho de los pulsos en 
funcion de la excitacion electrica hasta Rae!vl,5, y ese hecho ha sido confirmado 
experimentalmente para los valores de Rae empleados en las experiencias, ya que los valores 
obtenidos con el ajuste de las senales no han evidenciado dependencia de la excitacion 
electrica. 

En resumen de lo expresado acerca del estudio sobre los limites del paquete de ondas, 
puede decirse que existe una importante discrepancia entre las observaciones experimentales 
y las predicciones teoricas en lo que respecta al ancho de los paquetes de onda (de un orden 
de magnitud). Este hecho probablemente este asociado a que en un sistema real es imposible 
producir una perturbacion infinitesimal tipo pulso, por lo cual tal vez el paquete de ondas 
posea una extension inicial finita (es imposible obtener un movimiento del piezoelectrico 
representado por un pulso infinitesimal cuando el mismo es excitado con una senal electrica, y 
si asi fuera, es imposible que el chorro responda al mismo desarroUando una perturbacion 
infinitesimal tipo pulso sobre su superficie, debido al efecto de la viscosidad). La expansion 
del paquete de ondas previsto por la teoria, a medida que el mismo se aleja del orificio, 
concuerda cualitativamente con las observaciones experimentales, y este caso el desacuerdo 
mostrado posiblemente se vincule a las dificultades ocasionadas por el ruido, aunque no 
deben descartarse efectos viscosos y no lineales. Para las condiciones experimentales de este 
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trabajo, no se han observado cambios en la extension de los pulsos en funcion de la excitacion 
electrica, si bien los mismos serian, en teoria, de aproximadamente un 10% para la velocidad 
mas baja empleada. 

4.5 - CONCLUSIONES 

Las limitaciones al estudio experimental, propias del montaje utilizado algunas, y 
asociadas con la dinamica del escurrimiento otras, han sido analizas en primer lugar. Es a 
partir de tal analisis que se ha justificado la eleccion de algunos de los parametros 
experimentales. 

Con posterioridad a ello, los procedimientos que se han seguido en el tratamiento de las 
seiiales, tendientes a la determinacion de las tasas de crecimiento maximas y del ancho de los 
pulsos, han sido detallados. 

Finalmente, se han incluido los resultados obtenidos en lo que respecta a la forma que 
muestran las senales caracteristicas, la extension de los paquetes de ondas y las tasas de 
crecimiento maximas. Estos resultados han sido analizados y comparados con las previsiones 
teoricas, y las conclusiones que de ello pueden extraerse, se resumen como sigue: 

En lo concemiente a la extension del paquete de ondas, los resultados muestran que 
existe una importante discrepancia entre las predicciones teoricas previamente desarroUadas y 
las experiencias. Esto puede estar causado por una perturbacion que en su inicio tiene una 
extension finita, lo cual constituye un alejamiento de la hipotesis de una perturbacion 
infinitesimal tipo pulso. No obstante ello, algunas caracteristicas de la forma que adopta la 
superficie del chorro perturbado, y la expansion que experimenta el paquete de ondas a 
medida que se aleja del orificio, son correctamente descriptas por la teoria, en forma 
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cualitativa. En las condiciones para las cuales las experiencias han sido Uevadas a cabo, no se 
han observado cambios en la extension de los pulsos cuando la excitacion electrica es forzada. 
Cuando este fenomeno estaba previsto teoricamente, eran esperados cambios menores, del 
orden del 10%. De producirse realmente este cambio, es posible que haya quedado 
comprendido en el ruido de la serial. 

Respecto del efecto del campo electrico sobre la tasa de crecimiento maxima, la teoria 
concuerda con la tendencia mostrada en los experimentos (salvo cuando el chorro es afectado 
por efectos gravitatorios). Valores menores que los teoricos han sido obtenidos 
experimentalmente, pero el comportamiento monotonico con el campo electrico ha sido 
observado, aunque en menor grado para excitaciones electricas fuertes. 

En cuanto al comportamiento inestable del chorro en funcion de la velocidad del 
escurrimiento, la teoria se ajusta cualitativamente a las observaciones experimentales. El 
chorro se muestra mas estable a medida que la velocidad va en aumento. 

Las discrepancias encontradas indican que el modelo deberia ser refinado considerando 
la viscosidad del liquido, los efectos no lineales, y otro tipo de perturbacion a la superficie del 
mismo, diferente a un pulso infinitesimal (una entrada tipo rampa de elevada pendiente que se 
interrumpa para una amplitud muy baja, se constituye en una probable opcion para ello). Un 
analisis no lineal permitiria ademas predecir la longitud de ruptura y las caracteristicas de las 
gotas formadas. 
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ANEXO 4.1 - GRAFICOS PARA LA EVALUACION DE TASAS DE CRECIMIENTO 
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El objetivo planteado para esta tesis ha sido alcanzado: un estudio experimental de la 
propagacion de perturbaciones pulsantes en un chorro Hquido sometido a campos electricos 
pudo ser Uevado a cabo. 

En este informe, se ha efectuado un tratamiento basico de la dinamica de los chorros 
Hquidos y de sus aplicaciones, en el capitulo 1 . Las variables que intervienen en el fenomeno 
de la inestabilidad de los mismos, junto con la forma en la cual lo hacen, han sido puestas en 
evidencia. Las diferencias y similitudes que existen entre los chorros liquidos no excitados y 
los que si lo son, en cuanto a los modos de la superficie perturbada, regimenes, y tipos de 
ruptura, fueron tambien tratadas en esa oportunidad, acompafiadas con la explicacion de las 
vias usuales de producir la excitacion forzada. 

En el segundo capitulo el tema del estudio de inestabilidades ha sido abordado con una 
mayor profundidad. Las formas de Uevar a cabo estos estudios, y las clasificaciones corrientes 
que se emplean para los distintos tipos de inestabilidades han sido tratadas en forma 
introductoria. Un trabajo teorico desarroUado con anterioridad a este se presento en esa 
ocasion, en forma resumida, dando un marco de referenda para la realizacion del estudio 
experimental. 

Un dispositivo experimental, de una precision relativamente alta, ha sido concebido y 
construido para observar el fenomeno en estudio. La presentacion de su principio de 
funcionamiento y de sus caracteristicas constructivas ha sido efectuada en el capitulo 3. 

Finalmente, en el capitulo 4 se ha dado lugar a la presentacion de los resultados 

experimentales, junto con su discusion y comparacion frente a las predicciones teoricas. 

Como resultado de lo expuesto a lo largo del mismo, surge que el modelo teorico que 
fuera presentado acerca del fenomeno que aqui se trata, presenta desviaciones sobre todo 



150 



Conchisiones 



desde el punto de vista cuantitativo. Esas desviaciones podrian ser significativamente 
reducidas o tal vez anuladas si se tomase en consideracion la viscosidad del liquido, otro tipo 
de perturbacion inicial y una no linealidad en el tratamiento matematico de la superficie del 
chorro perturbado. El desarroUo de una teoria no lineal permitiria ademas predecir la longitud 
de ruptura y las caracteristicas de las gotas que surgen de la ruptura. 

Desde un punto experimental, en estudios futuros podrian ser consideradas algunas 
mejoras. Una mejora que surge de forma inmediata es la modificacion del sistema de medida 
del desplazamiento axial del chorro, ya que el uso de uno mas preciso reduciria el efecto de 
una de las fuentes de incertidumbre en los resultados. El sistema de presurizacion del 
recipiente de liquido podria ser revisado, ya sea empleando un sistema de nivel constante, o 
bien un sistema de control de la presion del gas inerte mas preciso que el utilizado. La 
formacion del chorro mediante un capilar dispuesto en el extremo de una camara, a la cual se 
pueda adosar el piezoelectrico (en este caso no perforado), podria ser considerada. De esa 
manera la medicion en posiciones mas cercanas al orificio podria ser practicada, aunque 
algunos efectos no deseados se presentarian (oscilaciones del capilar en voladizo y reflexiones 
de las ondas de presion dentro de la camara). Seria interesante por sobre todo, estudiar el 
fenomeno para valores mas bajos de We. Para ello, el radio del chorro deberia ser aumentado 
si se mantienen las mismas velocidades, pero los efectos gravitatorios se harian mas 
pronunciados y la limitacion por descarga corona se veria potenciada reduciendo el estudio a 
valores mas bajos de Rae. 

Respecto de las aplicaciones, sigue abierta la posibilidad de gobemar a voluntad la 
transicion entre el regimen de un chorro establecido en su configuracion inicial y el regimen 
de goteo, mediante la excitacion electrica del escurrimiento. No se han obtenido resultados 
categoricos respecto de la modificacion de los limites del paquete de ondas, en ese sentido. 
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por limitaciones propias del estudio experimental. En el unico caso analizado 
experimentalmente, para el cual se esperaban cambios (aunque pequeiios) en el ancho de los 
pulsos, ninguna modificacion ha sido observada tal vez debido al ruido en las senales 
adquiridas. No obstante, la concordancia cualitativa entre la teoria y las experiencias, que se 
ha mostrado para el resto de los parametros de la evolucion de estas inestabilidades, permite 
decir que un dispositivo que posibilite el cambio entre inestabilidades absolutas y convectivas 
en un chorro liquido, por medio de la excitacion con campos electricos, tal vez pueda ser 
desarroUado. 

Otra aplicacion de los resultados de este trabajo, se abre en la posibilidad de efectuar la 
medicion de caudales en un chorro liquido, en forma continua y poco invasiva para el 
escurrimiento, lo cual puede ser de gran interes en situaciones para las cuales se necesita una 
gran precision y haya transitorios importantes. En efecto, podria idearse un sistema 
constituido en parte por dos sensores ubicados en diferentes posiciones axiales, que permitan 
determinar la amplitud maxima del paquete de ondas, a traves de ello la tasa de crecimiento y 
por ende, de manera indirecta, la velocidad del escurrimiento. 
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